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Executive Summary 

The fatigue crack growth resistance of aerospace alloys is generally reduced by 
concomitant cyclic plastic deformation and exposure to a wide range of aggressive 
environments at ambient and elevated temperatures. Notable examples are iron and nickel- 
based superalloys in high pressure gaseous hydrogen, high strength alloy steels in water 
vapor, as well as aluminum, titanium and ferrous alloys in aqueous chloride solutions. 
Despite substantial advances over the past 30 years, the fracture mechanics approach to 
damage tolerant fatigue life prediction has not been adequately developed to enable 
incorporation of environmental effects on fatigue in codes such as NASA FLAGRO. 

Research, conducted at the University of Virginia (UVa) over the past four years 
and under NASA grant sponsorship, has aimed to develop the data, methods and 
understanding necessary to incorporate descriptions of environment-enhanced fatigue crack 
propagation (FCP) into NASA FLAGRO. The first phase produced a literature review of 
monotonic load stress corrosion cracking and environmental FCP (or corrosion fatigue) of a 
variety of structural alloys and environments relevant to aerospace applications.' The 
second phase of this research included three projects, conducted by undergraduates in the 
Department of Mechanical and Aerospace Engineering at UVa to fulfill the thesis-project 
requirement of the BS degree. This report documents the results of this team research. 

One student, Edward Richey III, has continued this work in conjunction with a Master of 
Science degree program in Mechanical Engineering. He was responsible for integrating the 
results of each undergraduate project into this report and the accompanying software 
package, and is continuing experimental work on corrosion fatigue crack propagation in the 
Ti-6Al-4V/aqueous chloride system. 

The underlying objective of the second phase of the research was to establish the 
foundation to incorporate deleterious environmental effects on fatigue crack propagation 
laws into NASA FLAGRO. Known methods to interpolate and potentially extrapolate 


R.P. Gangloff and S.S. Kim, Environnie/it Enhanced Fatigue Crack Propagation in Metals, NASA 
Contractor Report 4301, NASA Scientific and Technical Inl'ormation Division, 1990. 
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environmental FCP data were emphasized, including linear superposition and empirical 
curve-fitting descriptions of fatigue crack growth rate (da/dN) versus applied or effective 
stress intensity range (AK). The third approach, mechanism-based predictions of 
environmental FCP, is beyond the scope of this project and is being examined elsewhere.^ 

The results of the undergraduate projects were incorporated into a single computer 
program.* This program contains elements that were extracted from the NASA FLAGRO 
source code, particularly library data on material fatigue and fracture properties, as well as 
the Forman equations for calculations of fatigue crack growth rate versus stress intensity 
relationships, including the effect of crack closure. The program was written in Fortran 
and is executed by entering MENU from the DOS prompt on an IBM-compatible personal 
computer. The user can then select one of the following three options that resulted from 
each undergraduate project. 

Digitization of Crack Growth Rate Data (Jonathan M. Pope) 

The objective of this task was to develop a method to digitize FCP kinetics data, 
generally presented in terms of extensive da/dN- A K pairs, to produce a file for subsequent 
linear superposition or curve-fitting analysis. The method that was developed is specific to 
the Numonics 2400 Digitablet and is comparable to commercially available software 
products such as Digimatic™ Experiments demonstrated that the errors introduced by 
the photocopying of literature data, and digitization, are small compared to those inherent in 
laboratory methods to characterize FCP in benign and aggressive environments. The 
digitizing procedure was employed to obtain fifteen crack growth rate data sets for several 
aerospace alloys in aggressive environments. 


R.P. GanglofF, Corrosion Fatigue Crack Propagation in Metals y in Environment Induced Cracking 
of Metals , R.P. Gangloff and M,B. Ives, ^s., NACE, Houston, TX, pp. 55-109, 1990. 

R.P. Gangloff, NASA-UVa Light Aerospace Alloy ami Structures Technology Program y UVa Report 
No. UVA/528266/MSE94/1I4, March, 1994. 

A copy of this program may be obtained from the NASA Grant monitor. Dr. R.S. Piascik, in the 
Mechanics of Materials Branch at the NASA-Langley Research Center. 

Digimatic^^ is available from Famous Engineer Brand Software located in Richmond, Virginia. 
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Interpolation Modeling of Environmental FCP Data (Edward Richey, III) 

The objective of this task was to develop curve-fitting procedures for several well- 
known relationships between da/dN and AK, and to employ algorithms to interpolate the 
effects of loading frequency, hold time and stress ratio on the constants in these crack 
growth rate equations. Equations include the simple Paris power-law, the Forman Equation 
(with and without closure as developed in NASA FLAGRO), the Hyperbolic Sine Equation, 
and the Sigmoidal Equation. The program allows the user to select the form of the crack 
growth rate equation, as well as which constants to fix and which to optimize based on an 
input da/dN-AK data set. Interpolation is accomplished based on a logarithmic function, 
however, additional research is required to improve this approach by incorporating known 
results on time-cycle-dependent environmental FCP. Once the equation constants are 
determined as a function of frequency, hold-time and stress ratio based on several sets of 
input data for a given material and environment system, da/dN versus AK relationships can 
be interpolated. 

This empirical approach was assessed for Ti-6A1-4V in aqueous NaCl at several 
loading frequencies. The Forman, Hyperbolic Sine and Sigmoidal equations each 
reasonably described non-power-law crack growth rate behavior. The program incorporates 
statistical parameters (the coefficient of determination, r^, and a confidence interval 
estimate for each constant) to evaluate the goodness of fit for each crack growth rate 
equation. On this basis, r^ was generally above 0.90 for each of these three equations 
applied to six da/dN-AK data sets for Ti-6AI-4V in air and NaCl at three stress ratios. The 
Forman Equation (with or without closure) yielded the highest r^ for three cases, the 
Hyp>erbolic Sine Equation yielded the highest r^ for two cases, and the Sigmoidal Equation 
yielded the highest r^ for three cases. The Forman Equation and the Sigmoidal Equation 
yielded identically high average r^ values for all data. The program contains the graphics 
capability to plot da/dN versus AK data with the sup)erimposed curve fits. 

Linear Superposition Modeling of Environmental FCP Data (Allen W. Wilson) 

The basic linear superposition model put forth by Wei and Landes was implemented 
in the computer program, with improved analysis features. This approach enables the user 
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to input material property data on fatigue crack propagation in a benign environment, as 
well as data on monotonic load stress corrosion cracking for the embrittling environment of 
interest, coupled with the time and stress intensity dependence of the loading function. 
Material property data can be either extracted from the NASA-FLAGRO library, or entered 
as specific da/dN-AK or da/dt-K pairs and fit to the various equations described above. 

The portion of the fatigue load-cycle where environment is damaging must be determined 
by the user. With this input, the linear superposition program provides an output file and 
graphical representation of environmental fatigue crack growth rate versus AK, as a 
function of loading frequency, hold time, waveform, and stress ratio. The accuracy of the 
program was confirmed by comparisons with published linear superposition calculations. 

The results of this research augment the capability of NASA-FLAGRO. Example 
calculations demonstrating each element of the program are presented in several appendices 
to this rep>ort. Additional research is, however, required. The linear superposition 
approach is only accurate and relevant to a limited class of aerospace alloys that are 
exceptionally prone to environmental crack growth’ under monotonic loading. Those alloys 
that are properly selected to resist such cracking are likely to be sensitive to environmental 
FCP, however, the basic linear superposition approach will substantially underpredict 
da/dN. The curve fitting approach provides reasonable interpolations of the constants in the 
crack growth rate equations, however, the form of the inteipolation functions must be 
explored and coupled with physical understanding of time-cycle-dependent environmental 
FCP. This approach is not well suited for extrapolations outside of the input data base that 
was used to define the equation constants. These undergraduate projects further confirm the 
necessity to couple engineering estimation with fundamental understanding of environmental 
fatigue damage mechanisms.^ 


R,P. Wei and R.P. Gangloff» Enviromneufally Assisted Crack Growth in Structural Alloys: 
Perspeaives and New Directions, in Fracture Mechanics: Perspectives and Directions, ASTM STP 
1020 . R.P. Wei and R.P. Gangloff, eds., ASTM. Philadelphia, PA, pp. 233-264 (1989). 
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Introduction 


Overview 

One major problem when designing or analyzing a structural component is 
determining the behavior under cyclic loads, which cause cracking by fatigue, particularly 
coupled with corrosive environments. NASA would like to obtain an effective method of 
estimating environmental effects on fatigue crack propagation (FCP) rates in metals for 
use in damage tolerant codes such as FLAGRO. The objective of this research is to 
accomplish this task through a group undergraduate project within the Materials Science 
and Engineering Department at the University of Virginia. The method that will be 
developed for NASA will estimate the rate at which a crack grows when a structural 
alloy is exposed to cyclic loads in a corrosive or embrittling environment. Data from 
previous tests will be analyzed. 


Fatigue Behavior of Metals 

One methodology that has been developed for fatigue analysis is linear elastic 
fracture mechanics (LEFM). LEFM assumes that all metal parts contain preexisting 
cracks of engineering size. The speed with which these cracks propagate depends on the 
type of load, operating conditions, and environment. Cracks grow in a slow, stable 
manner until a critical crack size is reached, the crack becomes unstable, and the part 
catastrophically fractures. 

Linear elastic fracture mechanics is primarily concerned with predicting conditions 


1 



under which unstable crack growth will occur. The theory revolves around a relation 
involving the nominal stress (a), crack length (a), and stress intensity factor (K). The 
stress intensity factor is defined as: 

K - a y/n a Y (D 

where Y is a correction factor which depends on 
crack geometiy. Failure of the part is assumed to 
occur when the stress intensity factor exceeds a 
critical value which marks the onset of fast 
fracture. When fast fracture occurs, cracks in the 
part unstably grow at the speed of sound in the 
material, causing catastrophic failure. The critical 
value is usually the fracture toughness of the 
material, Kc, or Kjc- 

A typical fatigue crack growth rate curve is 
shown in Figure 1 . The curve is divided into 
three regions. Region I is the near threshold region where a decrease in AK (AK = 

- Kmin) results in a sharp decrease in da/dN as the threshold value AK* is reached. 

Cyclic loads which produce AK less than AK* do not cause crack propagation. Region 
in is the near critical region where a small increase in AK produces a large increase in 
da/dN as K„,ax approaches Kc. Region II is a transition region between Regions I and 
in, and can often be characterized by a simple power-law relationship between da/dN 
and AK. 



Figure 1: Schematic of three regimes 
of fatigue crack growth *. 
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Environmental Effects on Crack Growth 


The rate at which fatigue crack 
growth occurs is strongly dependent on 
the environment in which the part is 
operated. An inert or benign 
environment does not affect fatigue crack 
growth. A corrosive or embrittling 
environment affects the fatigue crack 
growth. 

In an inert environment, the rate 
of fatigue crack growth, (da/dN)f,tig„e> is 
generally only a function of the stress 
intensity factor range, AK, and the stress 
ratio, R (R = K,„i„/K^). Load 
characteristics such as frequency, wave 
form, and hold time do not greatly affect 
the behavior. In a corrosive environment. 





Figure 2: Effect of frequency on crack 
growth in Ti-8Al-lMo-lV alloy '. 


the load characteristics do affect the rate of 
crack growth. Frequency becomes an important factor because at lower frequencies 
"more time is allowed for environmental attack during the fatigue process.'” This effect 
is illustrated in Figure 2. The lower the frequency and the longer the hold time, the 
more time the environment has to attack the metal during each load cycle, thus the 
faster the crack grows. Additionally, environmental fatigue crack propagation 


can 
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involve complex dependencies of da/dN on AK and R; this behavior is unique compared 
to that typical of inert environment crack growth. 

Environmental subcritical crack growth occurs in corrosive environments without 
necessarily involving cyclic loads. The stress intensity is held coilstant and a time-based 
crack growth rate for the metal in the corrosive environment is determined. This crack 
growth rate has the form (da/dt)e„vjro„n.ent. which is the change in the crack length, a, per 

unit time, t (da/dt = 0 in a vacuum and without crack tip creep damage). There is no 

relevance of stress ratio or frequency since there are no cyclic load patterns and the 
stress intensity, K, is constant. Time is the main concern. 

Stress corrosion describes the growth of cracks due to static loads on a metal in 
environments which enhance crack growth. The static load causes localized plasticity at 
the crack tip. A chemical reaction between the material and ite environment causes crack 
propagation when none would occur in the absence of the environment. Because of the 
chemical reaction at the crack tip, the term corrosion is used to describe this mode of 
crack growth, although the surrounding material may not show typical signs of corrosion^. 

Figure 3 is a graphical depiction of the power laws that dictate the static load 
environmental crack growth rate. 

( da/dt) ^ D Ry below K* ( 2 ) 

{da/dt) environment ^ B above K* ( 3 ) 

where B, D, y, z and K* are material constants. 

Experiments can give accurate estimates of the material property coefficients in 
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Figure 3: Static load 
environmental crack growth. 

K* is the critical stress intensity 
between the two power laws for 
time-based crack growth in a 
corrosive environment. 


the static load environmental crack growth rate laws. However, simple power laws do 
not apply when both fatigue crack growth and static load environmental crack growth 
occur simultaneously. This situation is known as environmental fatigue crack growth (or 
corrosion fatigue), and the rate of crack growth is represented by (da/dN)totai- Several 
factors are pertinent in describing the conditions of environmental crack growth. These 
factors are the stress ratio (R), the frequency of the applied load (f in Hertz), and the 
environment. Other test factors which have an effect on crack growth rate, but are not 
addressed in this paper, are environment chemistry, temperature and metallurgical 
variables. 

(da/dN)ioiai is not generally described by a simple power law. For example, Wei 
and Landes^ outlined a method for predicting (da/dN)„,^ based on superposition of the 
effects of the crack growth rates, (da/dN)fa,ig„e and (da/dt)e„virx.nmem . to yield the total crack 
growth rate after integration: 


log da/dt 
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(da/dW)totai = ^ (4) 

T 

where for example: 

= ii^i ’ + (^--Kmin) (l-COSO)t)/2 (5) 

with cj = 2 IT f (rad/sec) 

7 = 1/ f (seconds) 

K(t) may be any time dependent function for a given application. A graphical depiction 
of superposition is given in Figure 4. The integration term in Equation 4 is often termed 
the stress corrosion crack growth rate, (da/dN)s„^ corrosion- The total crack growth rate, 
(da/dN),o^ in Equation 4, is both time and cycle dependent. Gallagher and Wei^ note 
that crack growth is dominated by environmental effects at low frequencies and by 
mechanical fatigue at high frequencies. 

(da/dN)u,m may also be described by more complex empirical equations such as 
the Forman, Hyperbolic Sine (SINH), and Sigmoidal Equations. These equations can all 
be fit to sets of (da/dN)ioui data and the constants in the equations determined by least 
squares regression. The fitted constants in the equations can be related to the load 
characteristics such as frequency, hold time, and stress ratio using empirical relationships. 
Using these empirical relationships, the equation constants can be determined for any 
combination of loading frequency, hold time and stress ratio for a given material- 
environment system. Once these constants have been determined for the desired load 
characteristics, (da/dN)u,tai can be calculated for any value of AK. 
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Figure 4: Schematic diagram illustrating the Wei and Landes method of linear 
superposition analysis, (a) Stress intensity spectrum in fatigue (top left); (b) rate of 
crack growth under constant load in a corrosive environment (top right); (c) 
environmental contribution to crack growth in fatigue (bottom left); (d) predicted 
environmental fatigue crack growth rate (bottom right)^ 

Project Oreanization 

Two different mathematical schemes were used to develop a method for 
predicting fatigue crack growth rates in corrosive environments. Allen Wilson developed 
a FORTRAN program to determine the crack growth rates with a linear superposition 
model which combines static load environmental cracking time-rate 
((da/dt)environmeni) ^^ta with inert environment mechanical fatigue cycle-rate ((da/dN)fatigue ) 
data. Edward Richey modified the NASA FORTRAN program FLAGRO which 
contains curve fitting functions that relate crack growth rates to AK and the stress ratio, 
R. Though both programs attempt to generate the same output, different input data and 
different calculative methods are employed. The results may not be similar because of 
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known deficiencies in each approach, at least for some material-environment systems. 
Jonathan Pope developed a method for collecting, organizing, and preparing data for use 
in the testing and development of the two computer programs. 

The main benefit of this project is the reduction of testing. If fatigue crack 
growth rates are required for a specific set of conditions that have not been tested, time 
consuming laboratory experiments are needed. The computer programs can supplement 
these experiments. The results of the computer programs will then be applied to 
different machine components that NASA uses so that the lifetime of the component can 
be accurately predicted. Expensive and exhaustive tests will not be necessary in 
obtaining a reasonable estimate of the component life. 

The following project report is divided into three sections. The first section 
discusses the digitization process for preparing data files. The second section discusses 
the interpolative computer model, while the third section discusses the linear 
superposition model. 



Section I: Data Collection and Preparation 


1. 1: Data Collection 

The two computer models for environmental FCP rates require different types of 
input data. The superposition model combines (da/dt)c„vironment and (da/dN)fetig„e to arrive 
at (da/dN)y^i. The interpolative model starts with several variations of (da/dN)y,jj, data 
and uses power law approximations to account for the differences in the variables that 
affect fatigue crack growth in a corrosive environment. The interpolative model can also 
describe (da/dt)e„yironnient and (da/dN)fajig„j data with equations. Therefore, all three types 
of crack growth rate data are needed to create a database for the two computer 
programs. These data will come from laboratory tests during the last twenty-five years. 
Each graph that is collected often has more than one set of useful data. For a specific 
alloy, an experimenter has varied either the environment, frequency of the load, or the 
stress ratio. 

Crack growth rate data are typically in the form of graphs with individual points 
specified and trend lines through the points, but only the data points are the results of 
laboratory testing. Unfortunately, the common practice is to only publish the graphs with 
the data points, while the precise numerical values of these data are required for the 
computer models. These values are often not retained by the experimenters, so 
searching for values from the authors of the tests would be a hopeless endeavor. As 
many as 200 test points may be represented on one graph, which makes it impractical to 
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handle all the numbers that would need to be published along with the graph. The 
intended outcome of collecting and digitizing a number of graphs is to provide a varied 
and complete set of data to be used in the computer programs. This was accomplished 
by recording copies of typical graphs that showed variations in the controlling 
characteristics on which the modeling programs would focus. The graphs were collected 
for specific aerospace alloys, environments, and variables. 

The database that is compiled is not a complete tool for crack growth rate 
prediction, but rather provides a realistic test of the computer models. For this reason, a 
large number of alloys are not examined. A more significant aspect of the database is 
whether or not it encompasses a wide array of variables. These variables include the 
stress ratio, frequency, yield strength, and environment. The alloys should be structural 
materials, fully characterized by fracture mechanics methods. Ti-6A1-4V, a common 
titanium alloy, and 4340 steel each fit these two criteria. Data for 7075 aluminum have 
also been recorded. 

Data from four different environments were collected. Air with 40% relative 
humidity, distilled water, a 3.5% sodium chloride solution, and argon make up the 
environments. Argon is used as an equivalent to vacuum for (da/dN) fa,jp„ data. These 
four environments are represented in both the Ti-6A1-4V and 4340 steel databases. 

Stress ratio variations are also included. 

One factor which affects fatigue crack growth in corrosive enviro nm ents is the 
frequency of the applied load. More so than any other factor, frequency is altered to 
generate different growth rate laws for a material-environment system. In other words 
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the metal, environment, and stress ratio are constant throughout a set of tests, but the 
frequency is changed. At lower frequencies, larger crack growth rates per cycle are 
experienced than for high frequencies. This occurs because the crack growth rate is 
measured in length per cycle; at a lower frequency there is a longer period of time that 
the metal is exposed to the corrosive environment for one cycle. (da/dt)e„vironmem plays a 
more dominant role in crack growth for this situation. 

In the 4340 steel database, yield strength is also important. 4340 steel has a large 
range of possible yield strengths controlled by tempering. A factor of 100 can separate 
the (da/dt)e„vironment ^ates for 4340 steels of different yield strengths. Therefore, careful 
attention has been paid to the yield strength of the specimens used in the various 
experiments. All graphs in the database are for 4340 steel with a yield strength of 1300 
to 1700 MPa. 

The graphs that have been collected are from a wide range of sources. Most of 
the data originated in papers in various journals. The graphs have since been duplicated 
and reprinted in other texts, compendia, and papers. In the data search, several sources 
presented graphs of the same data. Sometimes the overall appearance of the graph was 
changed while others looked exactly alike. Searches for specific data were often directed 
at journal articles, but the most productive manner of acquiring considerable amounts of 
data was to use compendiums. The Atlas of Fatigue Curves^ was very helpful, as were 
the Aerospace Structural Metals Handbook *^ and the Structural Alloys Handbook ^. 

A complete list of the data file names and the characteristics that describe each 
file are given in Appendix A. These files have been sent through the digitization process 
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that is described in the following text. 


1.2: Digitizing the Graphs* 


The database consists of fourteen graphs that have the same basic form as that 
shown in Figure 5. This plot has a large number of data points for each curve. The 


different curves have different characteristics and each need an individual data file. 



Figure 5: Typical da/dN vs. AK graph ®. 


* A commercial program, DIGIMATIC, has recently been marketed for PC operation 
under Windows. This product presumably provides similar results to the approach 
described here but was not tested in this study. 
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The computer models need input in the form of X and Y coordinates for each da/dN 
and AK. One could try to read each point off the graph with a ruler or a similar 
technique, but that would be foolhardy. The results would be relatively inaccurate and 
the process would be tedious and time consuming. A better method is needed. 

The method of choice is to electronically read the data points. This can be done 
using the Numonics 2400 Digitablet. Its primary use is associated with map reading and 
architectural dimensioning, but it meets our needs. The digitablet has an 
electromagnetic mechanism with a hand controlled cursor and keypad.^ A 
microprocessor and display console are also part of the digitablet system. The digitablet 
has many graphical functions, but it is used here only to accept data points. The origin 
of the graph is set as a zero point to enter data. The user then clicks on each data point 
to send it to a data file. Unfortunately, absolute distances, which are correlated to 
graphical data values, are measured only in centimeters or inches. Another problem 
arises from the set up of the graphs. The scales of the graphs are almost always 
logarithmic, but the digitablet is not capable of handling log scales. The data that are 
accepted by the digitablet in the form (X,Y) are meaningless when they stand alone. 
Section 1.3 corrects these problems and creates accurate (da/dN)-AK or (da/dt)-K data 
files. 

A step by step procedure for using the digitablet is included in Appendix B. It 
covers several necessities such as properly aligning the axes, setting the origin, and 
selecting the reference points. A steady hand is required for good results, but the 
manual work that needs to be done in entering a set of j)oints is minimal. Setting up the 
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system to receive the data is more complicated. 

When using the digitablet, an editor must be initiated by an external source and 
the local source where data are to be retained. It is suggested that the UNIX system 
and the JOVE editor be used to receive data.’ The manual in Appendix B provides a 
set of guidelines. A user familiar with other editors or systems can work out an alternate 
method. The local source mentioned above is the computer terminal to the right of the 
digitablet. The external source can be any network computer. 

Assuming that UNIX is the system of choice, the first step is to get a UNIX 
account. Becoming familiar with the UNIX system can be accomplished with the help of 
Academic Computing Center documents. It is important to note that a UNIX 
account can be logged onto from more than one machine at the same time. This feature 
will be needed to receive data from the digitablet. 

The external source opens the UNIX session, and the local source connects to the 
session. This order of operation must be followed or the data file that is created will not 
be recorded. Next the JOVE editor is evoked only on the digitablet terminal." The 
data points can now be entered with the keypad. Digitablet values will appear on the 
screen with other messages. The strange format of this information should be ignored. 
The data will be in a neat format when the acquisition process is complete. 

After data entry is finished and the file has been saved, JOVE is terminated. The 
file now appears on the external source in a sensible form, a column of X values and a 
column of Y values. Once the file appears on the external source, the digitablet 
terminal can close its UNIX session. 
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Attention now turns to moving the data file to a place where it can be used. The 
FTP mode (File Transfer Protocol) can be entered while still in UNIX*^, and used to 
export the UNIX file to DOS. The A:, B:, or hard drive should be specified as the 
destination for the file. After the transfer, a data file consisting of two columns of 
numbers exists on a PC disk. The process of assigning real values to these data lies 
ahead. The manual in Appendix B states the specific steps needed to prepare for this 
transformation process. 

The end result of the digitization and transformation process is the values of the 
graph points. With the values of the points, it is possible to recreate the graph seen in 
Figure 5, as illustrated in Figure 6. 


1.3: Data File Transformations 

The transformation process is necessary because the Numonics 2400 Digitablet 
does not allow the user to set the scale of the axes. All Digitablet measurements are 
calibrated in inches or centimeters relative to an input origin, rather than the values on 
the axes of the graph that is being digitized. The Digitablet output needs to be 
transformed to the correct values whether the graph has linear or logarithmic scales. 

Originally, four FORTRAN programs were written to transform the raw data files 
into accurate, meaningful data files. The four programs were similar in setup and 
mathematical logic. The differences between the programs dealt with the orientations of 
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the axes. Different combinations of logarithmic and linear scales on the X and Y axes 
created the need for more than one program. 

An error analysis of the digitization and transformation process is presented in 
Section 1.5. One result of a preliminary error analysis was the exposure of inadequacies 
in the four FORTRAN programs. One executable program which prompts the user for 
information about the graph was developed to replace the other programs. Either 
logarithmic or linear scales on any or both axes can be handled by the program 
DIGITIZE.EXE. All the functions that the four FORTRAN programs were capable of 
performing are included in the single executable program. 

From the original plot, the user must input the following information for each 

axis: 

1. Logarithmic or linear scale. 

2. The graph value at the Digitablet origin. 

3. The graph value at an arbitrary reference point far from the origin. 

4. The absolute distance between the origin and the reference point as 

displayed on the Digitablet console. 

Two purposes are served by recording information about the graphs. First, the 
graph origin is related to the zero point of the Digitablet output. For linear scales this is 
a number that is added to the digitized numbers, but this would also be a multiplying 
factor if the scales were logarithmic. The numbers in the data file are given meaning by 
scaling the axes. This is done with two values. One number is the graph value at a 
reference point, and the second number is the distance to the reference point as 
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reported on the Digitablet console. For example, on a log axis, the following data might 
be recorded. 

Origin: Digitablet reads 0, graph reads 10.0 ksi/in 

Reference point: Digitablet reads 6.000 cm, graph reads 80.0 ksi/in 

If the digitablet value of 4.0 cm is entered for a point, x, the transformation is 

X = (10.0) (80.0/10.0)^'* ®^®°^ = a graphical value of 40.0ksi\/in 

A similar function is set up to calculate values for linear scales. For the same graph data 
as in the example above, 

X = 10.0 + (4.0/6. 0)(80.0 - 10.0) = a graphical value of 56.67 ksiv^in 

The digitization manual in Appendix B includes a description of all the 
characteristics of the graph that need to be recorded. The program source code in 
Appendix C, DIGITIZE.FOR, is easily modified, so the user may set up the 
transformation in the manner that will yield the most accurate output for different 
situations. Adaptations to the programs should be considered for applications other than 
crack growth rate data. For typical crack growth rate graphs, the program will run 
accurately as it exists. 

Section 1.4 will move step by step through the digitization process beginning with 
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a graph. Section 1.5 discusses the error that can be expected in the results. These errors 
come not just from the reference point readings, but from many sources such as 
photocopying of graphs, shifting of the graph during digitization, and other human errors. 


1.4: Digitization and Transformation Example 

To illustrate the digitization and transformation process, one of the graphs from 
the database is presented as an example. One set of data is selected from the graph, and 
the results from the digitizer and 
transformation program are shown. All 
four sets of digitized data will be re- 
plotted at the end of the example. The 
original literature graph is shown in 
Figure 7, and the specific series that will 
be the focus of the example is defined 
by the inverted triangles. The defining 
characteristics of the data are: 
stress ratio =0.1 
frequency = 1 Hertz 
environment = 0.6 M NaCl 
units are in/cycle and ksi>/m 

Figure 7: Crack growth data for Ti-6A1-4V*‘*. 
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For the programmers who will use the data from these graphs, a list of the data 
files and their characteristics was created. This list appears in Appendix A. The data 
file which corresponds to the points in this example is TT1D.DAT. The first T in this 
name stands for transformed, the second T for titanium, the 1 for the first titanium 
graph, and the D for the fourth data set in the graph. 

During the digitization process, the information that is needed for the 
transformation process was recorded. The information that was recorded was: 


Digitablet origin: 

X axis value = 8 ksi-v/in 



Y axis value = 10 ’ in/cycle 


Reference point: 

X axis value = 150 ksi-v/in 

digitizer console = 6.492 


Y axis value = 10 ’ in/cycle 

digitizer console = 10.375 


All console units are in centimeters. The basis for choosing centimeters was to have the 
units consistent through all the digitized graphs. The following data file is the output 
from the digitizer. X Y 


1.038 

3.333 

1.025 

3.420 

1.270 

3.455 

1.528 

3.428 

1.498 

4.135 

1.783 

5.467 

1.793 

6.155 

1.845 

6.515 

2.010 

6.628 
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X 

Y 

2.113 

6.645 

2.135 

6.765 

2.050 

6.800 

2.403 

6.898 

2.433 

6.943 

2.313 

7.082 

2.653 

7.153 

2.715 

7.262 

2.963 

7.662 

3.060 

7.733 

3.153 

7.855 

3.448 

8.125 

3.635 

8.340 

4.057 

8.645 

4.598 

9.735 

-3.170 

-2.843 


The last two numbers in this series are negative to signal the end of the data to 
the transformation program. This point is mentioned in the digitization manual. This 
set of numbers was saved in the data file RT1D.DAT where the R stands for raw as 
opposed to transformed. 

The data file RT1D.DAT is ready to be transformed. DIGITIZE.EXE is run with 
the information that was taken from the graph. Both axes have log scales. The only 
other input for the program is the filename for the transformed data file. After 
execution, the output file contains the finished set of data seen below. 
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AK 

£ksi\/in) 

(da/dN)tocai 

fin/cvcle) 

12.782 

4.038 E-06 

12.708 

4.448 E-06 

14.194 

4.624 E-06 

14.571 

4.488 E-06 

15.733 

9.835 E-06 

17.894 

4.312 E-05 

17.975 

9.253 E-05 

18.402 

1.380 E-04 

19.825 

1.564 E-04 

20.768 

1.594 E-04 

20.976 

1.821 E-04 

20.186 

1.893 E-04 

23.674 

2. no E-04 

23.997 

2.218 E-04 

22.731 

2.588 E-04 

26.502 

2.800 E-04 

27.255 

3.161 E-04 

30.484 

4.926 E-04 

31.849 

5.330 E-04 

33.214 

6.103 E-04 

38.638 

8.235 E-04 

41.289 

1.045 E-03 

49.955 

1.466 E-03 

63.776 

4.916 E-03 


Since the original values of the data are not known, there is no way of gauging the 
accuracy of the digitization and transformation processes. The data can be checked by 
comparing the points on the graph to these values, but it is difficult to read numbers off 
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the axes. Another method is to replot the transformed data and compare the shapes of 
the curves. Figure 8 is the graph of the digitized data points. 

The differences between the original and the second version are small; no 
noticeable discrepancies can be identified when specific points are compared. Section 1.5 
establishes a range to the expected errors associated with the data acquisition technique. 


1.5: Error Analysis 

While digitizing and transforming a graph from the database, it became apparent 
that there is no way to quantify the error since the original graph values are not known. 
In an effort to define the accuracy of the method, a test has been conducted. The results 
of the test will give a range of error for the entire process which can then be compared 
to the typical uncertainty that comes from laboratory experimentation. 

A set of (da/dN)jotai data for aluminum alloy 2024 in moist air was acquired in 
both graphical and numerical form.*^ The graph was digitized and transformed, and the 
numbers from the transformation process were compared with the values of the original 
data. Figure 9 shows the original graph and Figure 10 shows a graph of the digitized and 
transformed data. Table I.l compares the results for specific points at the beginning and 
end of the graph. The average errors are for all the points in the data file 
(approximately forty). 

The data supplied to the modelling programs were often photocopied from 
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da/dN (in/ 


Ti-6A1-4V 

R=0.1 


vd 



5 7 10 30 50 70 100 


A K (ksiVin) 

Figure 8: Digitized and transformed set of data from Figure 7. 
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1 3 5 7 10 30 

A K (IfPaVm) 


Figure 9: Plot of data obtained in digital form for aluminum alloy 2024 tested in 
moist air at a constant K„^. 
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da/dN (mm/cycle) 


Digitized and Transformed Data 



1 3 5 7 10 30 

A K (MPaVm) 


Figure 10: Digitized and transformed set of data from Figure 9. 
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Table 1. 1: Error analysis for original graph (Figure 9) 





X original 
MPaVm 

X measured 
MPa/m 

3.26197 

3.26090 

3.29374 

3.28984 

3.32718 

3.32704 

3.39286 

3.38745 

3.47126 

3.47286 

3.58530 

3.59375 

3.71816 

3.73321 

9.11531 

9.13915 

9.38093 

9.41333 

9.70218 

9.76138 

9.91657 

9.97539 

10.2474 

10.2879 

10.5504 

10.5354 

10.8113 

10.8707 

11.0926 

11.1483 








0.160% 


0.046 % 


0.236% 


0.405 % 


0.262 % 


.345% 


0.610% 


0.593 % 


0.395 % 


.142% 


0.550 % 


0.502 % 


Y original 
mm/cycle 


9.77C-8 


1.19e-7 


1.31C-7 


1.71e-7 


.65e-7 


.96e-7 


8.41e-7 


Y measured % error 
mm/cycle 


6.54e-5 


7.83e-5 


8.50e-4 


.21e-5 


9.84e-5 


1.06e-4 


1.17e-4 


9.08e-8 


1.18e-7 


1.30e-7 


1.73e-7 


2.69e-7 


4.97e-7 


8.54e-7 


6.25e-5 


6.68e-5 


7.95e-5 


8.63e-5 


9.32e-5 


9.87e-5 


1.06e-4 


1.19e-4 


11 


IS 


IS 


.273 % 


. 866 % 


1.03% 


.785% 


1.54% 


.172% 


1.58% 


1.91 % 


2.19% 


1.53 % 


1.54% 


1.25% 


.257% 


.252% 


.986% 


IS 
























































































































different sources for use with the digitizer. For the most part, the source that the graph 
came from was not an original; this source probably copied the graph from another 
paper. There is no way of knowing how many times the graph has been reproduced, but 
the fact that the data are not first hand must be taken into account. To assess the effect 
of the copying of data, the test was run again with one variation. The original graph 
(Figure 9) was photocopied at 75% original size and then photocopied again at 133%. 
The graph is now approximately the original size. Figure 1 1 shows the re-plotted results 
of the photocopied data. The error analysis for the individual data points is presented in 
Table 1.2. 

The results in Tables I.l and 1.2 are similar when relating the errors on the X 
(AK) axis to errors on the Y (da/dN) axis. In both cases, the X axis errors are 
negligible. The largest error in either graph is less than 1% for the X value of the 
digitized point. The Y average error in Table I.l is 1.23%, and it increases to almost 
twice that value in Table 1.2. This suggests that Y axis errors are larger than X axis 
errors, and that photocopying magnifies the errors that result. 

Two questions must be addressed with respect to the errors presented in Tables 
I.l and 1.2. First, why are there errors? Second, why are there larger errors for the Y 
axis values than the X axis values? 

The answer to die first question goes back to the manual functions that are 
performed when digitizing the graphs. The Numonics 2400 Digitablet has a fixed 
coordinate system; when a graph is digitized, the axes of the graph must be aligned with 
the internal axes of the Digitablet. This requires patience since small movements create 
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da/dN (mm/cycle) 



1 3 5 7 10 30 

A K (llPaVm) 


Figure 11: Plot of data that were photocopied numerous times. 
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Table 1.2: Error analysis for photocopied graph (Figure 11) 


Y original 
mm/cycle 

Y measured 
mm/cycle 

9.77e-8 

9.30e-8 

1.19e-7 

1.22e-7 

1.31e-7 

1.35C-7 

1.71e-7 

1.77e-7 

2.65e-7 

2.74e-7 

4.96e-7 

5.14e-7 

8.41e-7 

8.69e-7 

6.13C-5 

6.29e-5 

6.54e-5 

6.71e-5 

7.83e-5 

7.94e-5 

8.50e^ 

8.58e-5 

9.21e-5 

9.28e-5 

9.84e-5 

9.83e-5 

1.06e^ 

1.05e-4 ( 

1.17e-4 

1.16e-4 ( 


|i^ 


X original 
MPaVin 


3.26197 


.29374 


3.32718 


3.39286 


3.47126 


3.58530 


3.71816 


9.11531 


9.38093 


9.70218 


9.91657 


10.2474 


10.5504 


10.8113 


11.0926 






X error 




Pa/m 


3.25497 


3.29016 


3.33587 


.39618 


3.47622 


.58054 


3.71651 


.06185 


.37736 


.70212 


9.96108 


10.2471 


10.5244 


10.7782 


11.0858 



.214% 


.109% 


.261 % 


.098% 


0.143% 


.132% 


0.044% 


.590% 


0.038 % 


.001 % 


0.447 % 


0 . 002 % 


0.247 % 


.306 % 


0.062 % 



maximum 

0.615% 

minimum 

0.001 % 

average 

0.240 % 


Y error 




.081 % 


.178% 


.700% 


.382% 


.457 % 


3.586 % 


3.280 % 


2.610% 


.654% 


1.384% 


.926 

% 

.785 

% 

.062 

% 

.981 

% 

.617 

% 
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noticeable errors in the digitized points. The user must continually adjust the graph until 
each axis strays no more than 0.01 to 0.02 centimeters in the perpendicular direction 
from the origin to the endpoint. 

Another manual error occurs when selecting points with the digitablet cursor. 
Points may be close together, or they may have odd shapes that make it difficult to find 
the center. It is easier to digitize a ’ + ’or an ’x’than it is to digitize a ’v’or a 
Finally, there is a tendency to move the cursor (mouse) when the enter button is pushed. 
Therefore, both axis alignment and the graph point marking suffer from human 
incapabilities. 

The next link in the chain is the transformation process. The mathematics of the 
program DIGITIZE.EXE are sound, but the input may be imprecise. A reference point 
on each axis is chosen, and the distance from the Digitablet origin to the reference point 
is recorded in order to scale the Digitablet output. Originally, the reference points were 
one cycle higher than the origin on logarithmic scales. This method proved to be 
ineffective when there were five full cycles present on the Y axis (da/dN data). Typical 
errors were 4 to 5%, but when the reference point was chosen to be as far from the 
origin as possible, errors in the Y axis values decreased to about 1 to 2 %. Improvements 
were noted for the X axis as well. Errors of 0.6 to 0.8% dropped to between 0.1 and 
0.3% 

There is no exact answer as to why the Y axis errors are larger than the X axis 
errors, but several possible explanations exist. The user is required to line up the axes 
with the Digitablet coordinate system. After photocopying a graph, the axes are often no 
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longer perpendicular, making it impossible to have both the X and Y axes correctly 
aligned. Often one axis is either consciously or subconsciously chosen to be more 
correctly aligned. Horizontal markings on the Digitablet surface coax the user to be 
more careful with the X axis than the Y axis. 

Another possible cause for greater accuracy in X values than in Y values is the 
structure of the data. The stress intensity on the X axis rarely varies more than one 
cycle on a logarithmic plot while the crack growth rate data on the Y axis may span two 
to six cycles. This wide range of values may degrade accuracy, since small errors will be 
magnifted by the large scale. This theory and the axis alignment theory were tested, 
both generating inconclusive results. 

The X and Y axes were switched in several tests to determine if axis alignment 
produced a noticeable difference in digitization errors. A sample of data was then 
retrieved from a file and plotted with just one cycle on the Y axis. The results suggest 
that the configuration of the da/dN data is more significant than axis alignment. The 
results of these tests along with operator bias tests are presented in Table 1.3. More 
photocopying examples are also included. 

All three operators achieved similar degrees of accuracy in digitizing the graph 
displayed in Figure 7. There was little variation in the average error for both the X and 
Y axes. Most users should be able to digitize graphs with accuracy comparable to the 
data presented in this table. 

The photocopying of the graphs presents another problem, distortion. The axes of 
a plot may not be perpendicular due to inadequacies in the copier, presenting an error 
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Table 1.3: Summary of test results 


Type of Graph 

Description 

avg. error 
in X value 

avg. error 
in Y value 


Log-log scale 

laser printer quality** 

0.209 % 

0.709 % 

- 


laser printer quality 

0.325 % 

1.225 % 



X and Y axes reversed 

2.13 % 

0.235 % 



(AK on Y and da/dN on X) 
data in one cycle on Y axis 

0.184 % 

0.294 % 



and one cycle on X axis 
photocopied twice 
(133 % then 75 %) 

0.240% 

1.97 % 


Semi-log 

photocopied once (1(X) %) 

0.128% 

1.087 % 



photocopied once (1(X) %) 

0.100% 

1.412% 



X and Y axes reversed 

0.290 % 

0.291 % 



one cycle on Y axis 

0.516 % 

0.418 % 


User bias test 

same graph for all three 




Log-log scale 

user 1 

0.325 % 

1.23 % 



user 2 

0.272 % 

1.17 % 



user 3 

0.541 % 

1.30 % 



source that can not be controlled. When this occurs, it is impossible to align both axes 
correctly at the beginning of the digitization process. The user must try to minimize 
these unavoidable errors on each axis. 


Plot went directly from graphics package to Hewlett Packard Laser Jet. 
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Based on the tests in Table 1.3, a conservative estimate of the errors involved with 
the data files is ±1% on the X axis and ±3% on the Y axis. The standard range of 
variation during laboratory tests can be much larger than these error estimates. A 
difference of 50% in da/dN for two specimens tested under the same conditions would 
not be unheard of. A 40% uncertainty band is presented for the test data in Figure 12. 
This estimate was created on the basis of other scatter bai^s presented for crack growth 
rate tests.’ All digitized data from Figures 10 and 11 easily fit in this band. Since the 
experimental data have a large uncertainty associated with fatigue crack growth rates in 
corrosive environments, the data that were supplied to the programmers should be 
sufficiently accurate to conduct realistic tests of the computer models. 


1.6: Conclusions for Data Collection and Preparation 

One purpose of this project was to create a database of crack growth rate data for 
use in two computer modelling programs. In creating this database, several techniques 
were developed to convert a graph with data points into an accurate file of digital data. 
The results may be useful to anyone interested in the process of digitization and 
transformation, as well as the programmers who will test their work here. 


1. A database of forty-six crack growth rate data files was creatol. The 
database includes static load environmental crack growth rate data 
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da/dN (mm/cycle) 


Scalterband and Original Data 



1 2 3 4 5 6 7 8610 20 30 


AK (llPaVm) 

Figure 12: 40% uncertainty band for test data. 
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((da/dt)j.„vironn^nt-K), incit environment fatigue crack growth rate data 
((da/dN)fjtigue-AK), and crack growth rate data for fatigue in a coirosive 
environment ((da/dN)„ia,-AK). Three alloys are included in the database 
which focuses on variations in environment, loading frequency, and stress 
ratio. 

2. A method for digitizing the points on a literature graph and transforming 
them to correct da/dN versus AK values was developed. Using the 
Numonics 2400 Digitablet, a data file is created by electronically reading 
graph points. The file is transferred to DOS and processed with a 
transformation program called DIGITIZE.EXE. The program converts the 
data to numerical values referenced to the axes of the original graph, for 
either linear or logarithmic scales. 

3. An error analysis was performed on the complete process of photocopying 
a graph, enlargement, digitization, and transformation. Photocopying 
increased errors, and human errors are unavoidable during digitization, but 
these effects were small. All errors are negligible when compared to the 
typical uncertainty associated with laboratory testing. The largest average 
error in any test was 2.1% while the smallest was 0.10%. The accuracy of 
the method is acceptable for supplying data for the computer models. 



Section II: Interpolative Model 


n.l: Crack Growth Rate Equations 

The interpolative computer model has the capability to fit a set of (da/dN)^,^ 
versus AK data to either the Forman Equation (with or without closure), the Hyperbolic 
Sine Equation (SINK), the Sigmoidal Equation, or the Paris Equation. The program 
allows the user to select the form of the crack growth rate equation, as well as which 
equation constants are to be optimized and which are to be fixed. The applicability of 
these equations to time dependent environmental FCP is assessed for Ti-6A1-4V in an 
aqueous chloride environment. 


Forman Equation (closure included) 

The Forman Equation including closure, developed by R.G. Forman, is discussed 
in detail in the NASA FLAGRO (Version 2.0)‘* user’s manual, and is of the general 
form: 


dN 


{ AlC V 

C AJC” 1--^^ 

V / 

il-R)^ [l--. — ^ Y 

[ ( 1 - 1 ?) K^j 


(II. 1) 


where N is the number of load cycles, a the is crack length, R the is stress ratio, and AK 
is the stress intensity factor range. The coefficient C, and the exponents n, p, and q, are 
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empirical constants which depend on material and environment. The parameter is 
the critical value of the stress intensity factor for unstable fracture, while AK^, is the 
threshold value of the stress intensity factor range. 

The crack opening function, f, is used to calculate the effect of stress ratio on 
crack growth for a constant amplitude load. The function was defined by Newman*® from 
crack closure considerations as; 


f = 


K 


op 




= max (i? i?^0 

A^+A^R -2^R<0 


(II. 2) 


where: 


A- = (0.825 - 0.34a + 0.05a**) 


cos 


^max \ 


I2 

jj 


Ai = (0.415 - 0.071a) 


A2 = 1 - A„ - Ai - A3 


(II. 3a) 

(II. 3b) 
(II. 3c) 


A 3 = 2A„ + A 3 - 1 (II. 3d) 

and Kop is the opening stress intensity factor, below which the crack is closed. The 
parameter a describes the condition of plane stress or plane strain, and 8 ^/ 0 ^ is the 
ratio of the maximum applied stress to the flow stress, where the flow stress is the stress 
required to produce plastic deformation in a solid. The parameter a is a material 
constant with values ranging from 1 for plane stress to 3 for plane strain. S^/Oo varies 
from experiment to experiment and is not a material constant, however, the value is 
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taken as 0.3 for most materials which were incorporated in the NASA FLAGRO 
program. Certain materials do not exhibit a large effect of stress ratio on da/dN, and 
the effects of crack closure do not need to be considered. For these materials a is set to 
5.845 and S^/Og to 1.0. The parameter f is then equal to zero for R 0 and becomes 
R for 0 ^ R < 1 . For positive stress ratios, the Forman Equation with closure then 
reduces to: 


dN 


( A ^ 

c AJC" 1 - ^ 
V AK , 


1 - 


AK 


(1 - R) Kr 


(II. 4) 


For negative stress ratios, AK is replaced by K,^, since only the positive part of a 
loading wave contributes to crack growth. Figure 13 shows a plot of f as a function of R 
for the cases represented by Equations II. 1 and II.4. 


Orisinal Forman Equation 

The original form of the Forman Equation, without an explicit treatment of crack 
closure (taken from NASA FLAGRO Version August 1986)”, is: 



(11.5) 


The coefficient C, and the exponents m, n, p, and q are all empirical constants. C, m, n, 
p, and q will not have the same values as in Equation II. 1 for a given set of data. 
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Figure 13; Crack closure function, f, as a function of stress 
ratio. From NASA FLAGRO User's Manual*®. 


Paris Equation 

The Paris Equation is a power law equation of the form: 

-^ = C (II. 6) 

ON 

where C and n are empirical constants. This equation effectively describes the FCP 
behavior in Region II, but since growth rates in Regions I and III are asymptotic to 
characteristic stress intensity ranges, this equation may not accurately model broad range 
da/dN versus AK data. 
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SiNH Equation 

The Hyperbolic Sine (SINH) Equation was originally developed by Pratt & 
Whitney, and is discussed by Haritos*® and Van Stone*’. The basic equation is of the 
form; 

= C^sinhCCj [log(A/0 + C 3 ] ) + (II. 7) 

where AK = K^-K^in and da/dN is the rate of crack growth per unit cycle. The 
parameters C3 and C4 are the horizontal and vertical locations of the inflection point 
respectively. The parameter C, is a constant with a value of 0 . 5 . The slope at the 
inflection point is defined to be C1C2. The SINH equation is symmetric about the 
inflection point, there is no mathematical definition of Kc or AK*. The curvature is 
identical in both the near toughness and near threshold regimes*’. The SINH Equation 
is schematically shown in Figure 14 . 

The parameters C2 through C4 are 
material-environment properties. To 
determine the values of the coefficients in 
the SINH equation, nonlinear regression 
must be used. This is done using an 
algorithm from Kuester^®. The algorithm 
uses a procedure proposed by MarquardF* to 
determine the parameters in a non-linear 
equation using the method of least squares. 



Figure 14; Schematic drawing 
of a SINH curve. 
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The procedure requires the user to input initial estimates for the parameters Cj, C3, and 
C4. These estimates can be determined from a graph of the data. Upper and lower 
bounds are then selected for each coefficient. Since the algorithm uses a least-squares 
objective function, the best value for each parameter is found by minimizing the sum of 
the squared errors for each data set. The sum of the squared errors, SSE, is defined as; 

SSE = 'L{ 

where Yp„d is the predicted value of da/dN and Yaaugi is the actual value of da/dN at a 
given AK. The algorithm determines the parameters by iteratively solving the following 
partial differential equations: 

-^^^ = -2E[p] [qcoshfCgdogAiC+C,)) (logAiC+Cj) ] =0 (il.9a) 

C/C2 

= -2 23 [p] [qcoshfCjdogAJT + C 3 ))]C 2 = 0 (ii.9b) 


dSSE 

dC^ 


-2S [pi = 0 


{II. 9c) 


P = [y-qsinh (C 2 (logAii: + C 3 ) ) - CJ (II. 9d) 

where y refers to the actual value of da/dN. 


Sismoidal Equation 

The Sigmoidal Equation was originally developed by GE Aircraft Engines, and is 


discussed by Van Stone*® and Haritos'* . The equation is of the form: 




(II. 10) 


AKc = (1 - i?) ifc 


where AKc is the critical value of the stress intensity factor range, and AKjh is the 
threshold value of the stress intensity factor range. A schematic of a sigmoidal curve is 
shown in Figure 15. The curve fitting 
parameters B, P, Q, D, AK^, and AKc are 
material-environment properties. 

The Sigmoidal Equation has an 
advantage over the SINK Equation 
because of deEned values for both AK,h 
and Kc. The Sigmoidal Equation is not 
identical in the near threshold and near 
toughness domains. Also, the Sigmoidal 
Equation deEnes the curvature in the 
transitions from Regions I and II, and 
Regions II and III. The main reason for this increased accuracy is the fact that the 
SINK Equation only has four curve Etting parameters, while the Sigmoidal Equation has 
six’’. The Sigmoidal Equation can also be rewritten in linear form, and the equation 
parameters determined using the least squares methods discussed below. 

To determine the equation parameters for the Forman, Paris and Sigmoidal 
Equations, the equations are Erst rewritten in linear form using logarithms. Using AK 



log A K 


Fignre 15: Schematic drawing of 
a sigmoidal curve . 
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and da/dN data, a set of simultaneous linear equations is constructed involving da/dN, 
AK, and the equation parameters. These simultaneous equations are solved using 
Gaussian Elimination and matrix algebra to determine the equation parameters. For 
example, the Forman Equation without closure (Equation II. 5 ) is rewritten as 

y = + A2X2 + ^3^3 + A^X^ + A5X5 (II. 11) 

where : 

(11.12a) 

Xj^ = log (1 - i?) (II. 12b) 

X 2 = log (AK) (11.12c) 


X 2 = log 



AK 



Xs = 1.0 


(II.12d) 

(II.12e) 

(II.12f) 


and Ai - m, Aj = n, A3 = p. At = q, and A5 = log C. Similarly, the Forman Equation 
with closure, the Paris Equation and the Sigmoidal Equation can be rewritten in linear 
form. 

The Forman Equation (with and without closure), the Paris Equation, and the 
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Sigmoidal Equation can be written in matrix form as: 


U] = [y] (11.13) 

where the A matrix contains the equation parameters such as C, n, p, q, and m, the X 
matrix contains the values of Xj through X5, and the Y matrix contains the values of log 
(da/dN). This equation can be rewritten as: 

[A] = [X]-^ [y] (11.14) 

The equation parameters can then be determined by Gaussian Elimination. As stated 
earlier, the SINK Equation can not be written in linear form, thus Gaussian Elimination 
can not be used for the SINH Equation. Instead, a non-linear least-squares algorithm 
must be used. 


Equation Parameters 

The FCP laws used in the interpolation model all change with stress ratio, 
frequency, and hold time variations. The interpolative model can be used to empirically 
describe such changes within the data range, but will not likely provide accurate 
predictions outside of the bounds of the experimental values. As a first approximation, 
the equation parameters, such as C, n, p, and q are set equal to linear relationships 
involving assumed logarithmic functions of the stress ratio, frequency and hold time of 
the applied load. The form of the relationship depends on the number of data sets 
available to the program and the form of the crack growth rate equation. The number 
of loading variables (R, hold time, frequency, K„^, and that the equation 
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parameters are related to is equal to the number of data sets entered. If only one 
da/dN versus AK data set is entered the crack growth rate equation parameters for a 
given frequency are set equal to the following relationship: 

q = log 1 1^1 + ij (11.15) 

f = frequency of the applied load 

where rj is an equation parameter such as C or n, and aj is constant fit by the computer. 
The user has the option of deciding which parameters to fit and which to fix. The 
accuracy of the fitted equation increases as the number of the parameters fitted 
increases. The form of the relationship was selected so that at low frequencies (higher 
corrosion effects), the log value increases for a fixed frequency. If two data sets are 
entered, and the user has selected the Forman Equation as tte crack growth rate model, 
the relationship is: 

q = log 1 + ^2 (11.16) 

f = frequency of the applied load 

T = hold time of the applied load 

If die user selects the Paris Equation, the SINH Equation, or the Sigmoidal 
Equation and enters two data sets, the relationship is: 

q = log 1^1 + 1+32 log [— (11.17) 
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f = frequency of the applied load 


R = stress ratio 

if the stress ratio is held constant. If is held constant for the data set then the 
relationship is; 

T) - log 1^1 + l] + 32 logli?;^) ( 11 . 18 ) 

If ^in Is hoW constant, in Equation 11.18 is replaced by K,„in- Th® form of the 
relationship was selected so that a dependency on the stress ratio could be introduced 
into the Paris, SINK and Sigmoidal Equations. The Forman Equation already describes 
the effect of stress ratio based on a crack closure strategy; that is why a hold time effect 
and not a stress ratio effect was included in Equation II. 16. 

The effects of hold time are not accounted for in the Paris, SINK or Sigmoidal 
Equations if two data sets are used. To take into account the hold time for the Paris, 
Sigmoidal or SINK Equations, three data sets must be used. If three data sets are 
entered, the following relationship is used for all of the crack growth rate equations: 

n = a,log |lj + ij * a,log * a,log[:^] ( 11 . 19 ) 

f = frequency of the applied load 
T = hold time of the applied load 
R = stress ratio 
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if the stress ratio is constant. If is held constant the relationship is : 


n 



+ 


^2 + ^3 log 


T + 1 

0.1 


( 11 . 20 ) 


If K„i„ is held constant, then in Equation 11.20 is replaced by Kn^n- These equations 
allow the user to determine the equation parameters for a particular data set if the load 
characteristics are known. The assumed logarithmic functions of f, R, K„ax. and r 
must be tested with experimental data. These functions should only be used for 
interpolation, not prediction of the effects of variables outside the range of the curve-fit 
data. The logarithmic forms of the relationships were selected so that the relationships 
were defined for the allowable ranges of frequency, hold time, K,^, K,^„, and stress ratio. 
For example, the expression involving the hold time is always deflned unless t = -1.0, 
which is not physically permissible. For all r ^ 0 the expression does not tend toward 
infinity. 

To determine the values of aj, aj, and a,, the equation parameters are determined 
for each data set entered. A set of simultaneous equations is then developed for each 
parameter which relates the parameter values to the given load characteristics. Matrix 
algebra and Gaussian Elimination are then used to solve these equations and determine 
the coefficients. 

Using Equations II. 15 through 11.20, the program allows the user to interpolate 
da/dN data given the stress ratio (or or frequency, hold time, and AK. This is 
done by first selecting the form of the crack growth rate equation, entering AK and 
da/dN data, and determining the values of the coefficients a,, aj, and aj in Equations 
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n. 15 through 11.20 for selected equation parameters. Given the value of the frequency, 
hold time, and stress ratio (or or K„,in), appropriate values for the crack growth rate 
equation parameters can be determined for any set of load conditions for the same 
material-environment system. Given a value of AK, da/dN can then be calculated. 


n.2: Interpolative Computer Modeling Program 
Computer Program Development 

The first step in developing the computer modeling program was to obtain a copy 
of the FLAGRO source code from R.G. Forman at the NASA- Johnson Space Flight 
Center. Many of the subroutines in the interpolative model were taken directly from 
FLAGRO. The next step was to locate the sections of the program that were essential. 
This was necessary because the FLAGRO program is too large to run on the 
FORTRAN compilers available at the University of Virginia. The sections that were 
removed to be used in the final program were the subroutines which loaded the material 
properties from a library data file for use in the crack growth rate equations, and the 
subroutines which controlled the entering of AK and da/dN data. 

The subroutines in FLAGRO which controlled the linear regression analysis used 
to determine the equation parameters offered the user a choice between the Forman 
Equation with closure or the Paris Equation. Also, the user was only given the choice of 
determining C, n, or both C and n in Equations II. 1 and II. 6. The code was first 
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modified so that the Forman Equation without closure, the SINK Equation, and the 
Sigmoidal Equation could be used to model da/dN data, as well as the existing Forman 
Equation with closure and the Paris Equation. The user is thus allowed to choose from 
five crack growth rate equations. The program was then altered so that the user could 
optimize any combination of the equation parameters, including the parameters C, n, p, 
q, and m. Next, subroutines which used least squares regression and Gaussian 
elimination to solve for the equation parameters were added to the program. The 
subroutines used to perform the least squares analysis and Gaussian Elimination were 
based upon FORTRAN programs written by Scott^ and Ellis^^. Subroutines to perform 
non-linear least squares regression were also added“. 

Subroutines to plot the actual and predicted da/dN data were also added to the 
program. These subroutines were based upon the Watcom Graphical Kernel System 
(WGKS) package. The actual input data were plotted as points, while the fitted crack 
growth rate equation was plotted as a smooth curve. This allows the user to judge how 
well the curve fit equation correlates with the input data. Finally, the program was 
modified so that the parameters which were optimized were set equal to the linear 
relationships (Equations 11.15 through 11.20) involving the stress ratio (or or 
frequency and hold time of the applied load. 

Comparison to NASA FLAGRO 

The stand alone interpolative computer program developed to model fatigue 
behavior contains many of the useful features found in the FLAGRO program. As 
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mentioned earlier, the materials files used in the FLAGRO program are incorporated 
into the stand alone program. When the FLAGRO program is used, however, the user 
is required to enter the material properties manually when fitting data to a crack growth 
rate equation. The option to fit data to an equation does not utilize the material files 
that are available. The stand alone program does utilize these files when the user so 
desires, and allows the user to enter the material properties from a material file. The 
material properties that are necessary are the yield strength, Kjc, A^, S^/a„, AK„, 

and a used in Equations 11.2,11.3, 11.23, and 11.24. These material properties are used 
when calculating AK^, Kc, and the crack opening function f. The equation parameters C, 
n, p, and q in the Forman Equation with closure are also loaded. These parameters are 
only loaded so that the user can see the values for these constants determined by NASA. 

Another difference between the stand alone program and the NASA FLAGRO 
program is the option to hold or constant for a set of data. The NASA 
FLAGRO program assumes that the stress ratio, R, is constant for the data entered. 

This assumption is not always valid since experiments are conducted with (or K^jn) 
being constant. If the user wishes to hold or K„,i„ constant, the program calculates 
the stress ratio for every point in the data set using one of the following equations: 

R = {II. 21) 

•^ax 


R = 


^in + Air 


( 11 . 22 ) 


This array of stress ratios is then used in the least squares regression used to determine 
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the equation parameters, and is also used when calculating the crack closure function, f, 
when necessary. The values of Kc and AK^, are constant for a given material data set. 

The values for and AK^, calculated by the NASA FLAGRO program are 
primarily accurate for fatigue in moist air. Since the stand-alone program was developed 
to model fatigue behavior in embrittling environments, these values are not always 
acceptable. An option to optimize Kc and AK^ has been incorporated into the stand- 
alone program in order to address this problem. The user presently has the option of 
entering Kc and AK^, manually, calculating the values from the equations contained in 
the NASA FLAGRO program using the material properties discussed earlier, or 
optimizing the values based on input da/dN versus AK data. If the user decides to 
calculate the value using the NASA FLAGRO equations, the following relationships are 
used:‘® 

AJiTti, - AJ«:„ /-ij tan-i (1 - i?) (11.23) 





(11.24) 


where t is the thickness, AK^ is the fatigue threshold at R =0, A,; and Bk are material 
constants, and Uys is the yield strength. If the user decides to optimize the values for Kc 
or AK*, the program fits the data set to a cubic equation using least squares regression. 
The cubic equation contains an upper and lower asymptote, and these asymptotes 
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correspond to the value of Kc and AKt^. These two parameters are then used in 
subsequent interpolative modeling using the Forman or Sigmoidal Equations. 

Several statistical parameters were also added to the stand alone program. First, 
the coefficient of determination, r^, is calculated for each data set fit to a crack growth 
rate equation. The coefficient of determination is the percentage of variation in the 
dependent variable, da/dN, which is accounted for by the relationship between da/dN 
and AK represented by the fitted crack growth rate equation. Mathematically the 
coefficient of determination is calculated as: 


2 - 2 = 1 -, 




-Y 

actual ■*‘'pied* 


S(K 


-Y 

actual average ' 


( 11 . 25 ) 


where Yacuai is the actual value of da/dN and Yp^ is the equation predicted value of 
da/dN at a given AK. Yav^age is the average value of the actual da/dN values. The 
coefficient of determination will range between 0 and 1.0. The closer the value is to 1.0, 
the better the fit of the crack growth rate equation. 

Confidence interval estimates for each fitted equation parameter were also added 
to the program. A 100(1 -a) percent confidence interval for the equation parameter 6j is 
defined as: 


where t is the student’s t distribution with n-p degrees of freedom. The quantity is 
given by: 


o 


2 


SSE 

n-p 


( 11 . 27 ) 
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and C|j is the jjth element of the ([X]’[X]) ‘ matrix. The [X] matrix is defined as before, 
there are n data points in the data set, and p equation parameters are fitted.^'* 


n.3: Data Testing and Results for the Interpolative Model 

Data Testing 

The first criteria used to test the interpolative computer model was to benchmark 
the programs’ linear least squares regression subroutines against the FLAGRO program. 
To do this, data set TTIB was curve fit to the Forman Equation with closure using the 
stand alone program and the FLAGRO program. The data was taken from an article by 
Dawson and Pelloux''*. The value for Kc was entered as 69 ksi\/in, AK^ was 7.27 ksi\/in, 
a was 2.5, and S^/Og was 0.3 for each program. The parameters C and n were fit, while 
p was 0.25 and q was 0.75 for each program. Results are shown in Table ILL 
Table ILL Bench marking of linear regression algorithm* 


Parameter 

FLAGRO 

Stand Alone 

Difference 

% Difference 

C 

0.5759 X 10’ 

0.5673 X 10’ 

8.6 X 10‘® 

1.49 

n 

2.271 

2.271 

0.0 

0.0 


‘The units for the equation constants depend on the units of the fitted da/dN-AK 
data. The fitted constants can only be used with data that have units consistent with the 
data used to determine the constants. 
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The results were identical for the parameter n, and varied only slightly for parameter C. 
Thus, the linear least squares regression analysis contained in the stand alone program is 
acceptable. 

The algorithm used to determine the parameters in the Hyperbolic Sine Equation 
was also bench marked using results from Haritos‘* for R = 0.1,f = 1 Hz, and t = 0.0 
seconds. The values selected by Haritos'* and the interpolative model are shown in 
Table II. 2. The values selected by the interpolative model did vary from those selected 
by Haritos*®, however, the values determined by the interpolative model yielded an of 
0.996, while Haritos’ values yielded an r^ of 0.992. Results are shown in Figure 16. 


Table II. 2: Bench marking of non-linear regression algorithm. 


Parameter 

Haritos 

Interpolative Model 

% Difference 

Cz 

3.197 

3.896 

21.9 

C3 

-1.762 

-1.6167 

8.2 

C4 

-2.530 

-2.8080 

11.0 


The second criteria used to test the accuracy of the computer model was to 
determine how well the program fit fatigue data to the various crack growth rate 
equations. The data used were for titanium alloy Ti-6A1-4V in both an air environment 
and a 1.0% NaCl solution. The data were taken from a NASA-UVA Light Aerospace 
Alloy and Structures Technology Progress Report^*. The frequency was constant at 5 Hz, 
the hold time was constant at 0.0 seconds for all of the data, and a was entered as 2.5. 
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Figure 16: Comparison of curves fit by Haritos^® and the interpolative model. 

The first equation tested was the Forman Equation with closure. Two tests were 
conducted for each data set. In the first test, only C and n were fit, while p and q were 
entered manually. The parameter p was set equal to 0.25 and q was 0.75. These values 
were taken from the NASA FLAGRO materials file. The values calculated by the 
program are shown in Table II. 3. In the second test, all of the equation parameters were 
optimized. The values selected by the program are shown in Table II.4. In both tests, 
the values for AK^ were determined from AK„, and the value of Kc was held constant at 
3953 MPa\/mm for air and 4270 MPav/mm for 1.0% NaCl. The value of AK„ was taken 
from the FLAGRO material files as 121.6 MPav^mm for Ti-6A1-4V (MA). 




The same tests were then conducted for the Forman Equation without closure, 
holding the parameter m constant at 1. The results for fitting only C and n are shown in 
Table II. 5, and the results for fitting C,n,p, and q are shown in Table II. 6. Finally, the 
data were fit to the SINH Equation and the Sigmoidal Equation; results are shown in 
Tables II. 7 and II. 8. Results are shown graphically in Figures 17 through 22. 

Notice that when fitting C, n, p, and q in the Forman Equation (with and without 
closure), as well as when fitting B, P, Q, and D in the Sigmoidal Equation, the equation 
is not always accurate outside of the fitted data range. The equation often becomes 
sinusoidal in nature when all four parameters are fit. This means that near the 
threshold, as AK decreases, da/dN increases. Also, as AK approaches AKp, da/dN 
decreases as AK increases. This is not physically permissible, thus when all four 
parameters are fit, the parameter values often cannot be used for AK values outside of 
the fitted data range. This problem can be corrected by fitting a larger range of data to 
the equations. If data from all three regions of crack growth are fit to the equation, 
from AK^ to AK^, this problem will be eliminated. Unfortunately, the data used in this 
study did not encompass all three regions of crack growth. 

Recently, Newman" suggested setting p and q equal to 2.0 when values are not 
definable from wide-range da/dN-AK data. The data for titanium alloy T1-6A1-4V in 


Pages 54 though 68 have been revised at the request of Dr. Robert S. Piascik, 
Mechanics of Materials Branch, NASA Langley Research Center, in order to include a 
third analysis case where p = q = 2.0. These values for p and q were selected by Dr. 
James C. Newman, Mechanics of Materials Branch, NASA Langley Research Center. 
This case was not included in analysis of model capability to interpolate trends in fatigue 
behavior. 
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both the air and 1.0% NaCl environments were fit to the Forman Equation with and 
without crack closure. Kc was held constant at 3953 MPav/mm for the air environment, 
and at 4270 MPa\/mm for the 1.0% NaCl solution. was entered as 0.3, a was 2.5, 

and AK,h was calculated from AK^. AK^ was taken from the FLAGRO materials files as 
121.6 MPa\/nun for T1-6A1-4V (MA). The values determined for C and n are shown in 
Tables II. 8a and II. 8b. Results are shown graphically in Figures 22a through 22f. 

Notice in Figures 22a through 22f, that setting p and q equal to 2.0 alleviates the 
problem encountered when fitting all four equation parameters in the Forman Equation. 
When holding p and q constant at 2.0, the equation does not tend to become sinusoidal, 
and the crack growth rates predicted by the equation are physically permissable, even 
outside the date range. The instability problem encountered with the Sigmoidal 
Equation has yet to be addressed. 


Table II. 3: Forman Equation (Equation II. 1) with closure (C,n Fit) 


Parameter 

Environment 

R = 0.1 

R =0.4 

R = 0.7 

c 

Air 

0.6032 X 10'’ 

0.8089X 10 '* 

0.1210X 10'* 

n 

Air 

5.036 

4.200 

3.762 

c 

1.0% NaCl 

0.1727X 10'* 

0.3308X 10'* 

0.3554X 10'* 

n 

1.0% NaC! 

5.028 

4.532 

4.095 
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Table II. 4: Forman Equation (Equation II. 1) with closure (C,n,p,qFit) 


Parameter 

Environment 

R = 0.1 

II 

o 

R = 0.7 

c 

Air 

0.1510X 10” 

0.2404 X 10” 

0.2829X 10” 

n 

Air 

22.07 

8.020 

8.129 

p 

Air 

-14.71 

-1.395 

-9.469 

q 

Air 

-78.72 

-16.43 

-6.056 

c 

1.0% NaCl 

0.4368 X 10” 

0.2289X 10” 

0.3345 X 10” 

n 

1.0% NaCl 

12.33 

-10.21 

-9.031 

p 

1.0% NaCl 

-6.494 

8.012 

16.43 

q 

1.0% NaCl 

-33.84 

99.19 

41.21 


Table II. 5: Forman Equation (Equation II.5) without closure (C,n Fit) 


Parameter 

Environment 

R = 0.1 

R =0.4 

R = 0.7 

c 

Air 

0.1917X 10” 

0.6918X 10 ” 

0.1673X 10” 

n 

Air 

5.036 

4.200 

3.762 

c 

1.0% NaCl 

0.5500X 10” 

0.2766X 10” 

0.4921 X 10'“ 

n 

1.0% NaCl 

5.028 

4.532 

4.095 
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Table II. 6: Forman Equation (Equation II. 5) without closure (C,n,p,q Fit) 


Parameter 

Environment 

R = 0.1 

R =0.4 

R = 0.7 

c 

Air 

0.3942X 10“ 

0.1029X 10” 

0.1453X 10-^® 

n 

Air 

22.07 

8.020 

8.129 

p 

Air 

-14.71 

-1.395 

-9.469 

q 

Air 

-78.72 

-16.43 

-6.056 

c 

1.0% NaCl 

0.9722X 10“ 

0.7540X 10^ 

0.2978 X 10 ‘ 

n 

1.0% NaCl 

12.33 

-10.21 

-9.031 

p 

1.0% NaCl 

-6.494 

8.012 

16.43 

q 

1.0% NaCl 

-33.84 

99.19 

41.21 


Table II. 7: Hyperbolic Sine Equation (Equation II. 7), C, =0.5 


Parameter 

Environment 

R = 0.1 

R = 0.4 

R = 0.7 

Cj 

Air 

2.8560 

1.4889 

1.0491 

C3 

Air 

-3.3127 

-4.1999 

-4.9851 

C4 

Air 

-2.6468 

-1.5061 

-1.0000 

C2 

1.0% NaCl 

1.4749 

1.7548 

2.1795 

C3 

1.0% NaCl 

-4.3175 

-1.0000 

-1.3854 

C4 

1.0% NaCl 

-1.0000 

-7.6946 

-6.8909 
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Table II. 8: Sigmoidal Equation (Equation II. 10) 


Parameter 

Environment 

R = 0.1 

R =0.4 

R = 0.7 

B 

Air 

-118.60 

-38.07 

-32.34 

p 

Air 

48.71 

13.60 

13.60 

D 

Air 

64.59 

14.15 

6.150 

Q 

Air 

-16.23 

-2.097 

-8.944 

B 

1.0% NaCl 

-56.96 

87.05 

46.99 

P 

1.0% NaCl 

22.31 

-34.64 

-28.53 

D 

1.0% NaCl 

25.40 

-72.27 

-37.41 

Q 

1.0% NaCl 

-6.652 

8.918 

17.51 


Table II. 8a: Forman Equation (Equation II. 1) with closure for p = q = 2 


Parameter 

Environment 

R =0.1 

R =0.4 

R = 0.7 

C 

Air 

0.1642X 10'' 

0.1317X 10'^ 

0.1029X 10” 

n 

Air 

4.144 

3.374 

2.968 

c 

1.0% NaCl 

0.1272X 10” 

0.6862 X 10'^ 

0.2976X 10” 

n 

1.0% NaCl 

3.970 

3.247 

3.310 
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Table II. 8b: Forman Equation (Equation II. 5) without closure for p = q = 2.0 


Parameter 

Environment 

R =0.1 

R =0.4 

R = 0.7 

c 

Air 

0.5630X 10’’ 

0.6290X 10 ’’ 

0.3141X 10'^ 

n 

Air 

4.144 

3.374 

2.968 

c 

1.0% NaCl 

0.4544 X 10'^ 

0.3305 X 10'^ 

0.9106X 10'^ 

n 

1.0% NaCl 

3.970 

3.247 

3.310 


In order to quantitatively determine how well the crack growth rate equations fit 
the data sets, the coefficient of determination was calculated for each equation and data 
set. These values are shown in Tables II. 9 and II. 10 for the original two tests conducted. 
Results for the case recommended by Newman of p = q = 2 in the Forman Equation 
(with and without closure) are shown in Tables II. 10a and II. 10b. 

The third criteria used to test the computer program was to see how well the 
program could interpolate trends in da/dN behavior. To accomplish this task, three sets 
of data for titanium alloy T1-6A1-4V in 3.5% NaCl solution were used; conditions were 
constant except for the loading frequency, thus all of the equation parameters had to be 
expressed as a function of frequency. The data were taken from a study by Dawson and 
Pelloux’^. To test the model, data for 1 Hz and 10 Hz were fit to the Forman Equation 
(with and without closure), the SINH Equation, and the Sigmoidal Equation, and the 
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C) SINK Equation. 



100 300 500 700 

A K (MPoV mm) 


D) Sigmoidal Equation. 


Figure 17; Fit of fatigue data for Ti-6A1-4V in air with R = 0.1 at 5 Hz. 








A) Forman Equation (closure). B) Forman Equation (no closure). 



Q SINH Equation. 



D) Sigmoidal Equation. 


Figure 18: Fit of fatigue data for Ti-6A1-4V in air with R = 0.4 at 5 Hz. 
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AK (MPoVmm) ^ ^ (MPoVmm) 


A) Forman Equation (closure). B) Forman Equation (no closure). 



70 100 300 500 70 100 300 500 

AK (MPaVmm) AK (MPaVmm) 

C) SINK Equation. D) Sigmoidal Equation. 

Figure 19; Fit of fatigue data for Ti-6A1-4V in air with R = 0.7 at 5 Hz. 
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T1-6AI-4V 

1 .0^ NoCI Environment 
-3 _ R = 0.1 

; Frequency = 5 Hz 


TI-6AI-4V 

1 .0% NaCI Environment 
_ R = 0.1 

: Frequency = 5 Hz 


Predicted (C,n fit) 

Predicted (C,n,p,q fit) ' 

O Actual Dota 


Predicted (C,n fit) 

Predicted (C.n, p,q fit) 

O Actual Dato 


300 500 700 

AK (MPaV mm) 


300 500 700 

AK (MPoV mm) 


A) Forman Equation (closure). 


B) Forman Equation (no closure). 


Ti-6AI-4V 

‘ 1 .0% NaCt Environment 
-3 R = 0.1 

r Frequency = 5 Hz 


Ti-6AI-4V 

1.0% NaCI Environment 
_-3 R = 0.1 
’ ^ : Frequency = 5 Hz 


O Actual Date 


Sigmoidol Fit 

O Actual Dota 


300 500 700 


A K (MPaVmm) 


C) SINK Equation. 


300 500 700 


A K (MPaVmm) 


D) Sigmoidal Equation. 


Figure 20: Fit of fatigue data for Ti-6AI-4V in 1.0% NaCI with R = 0.1 at 5 hz. 
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A) Forman Equation (closure). 


B) Forman Equation (no closure). 



C) SINK Equation. 



D) Sigmoidal Equation. 


Figure 21: Fit of fatigue data for Ti-6A1-4V in 1.0% NaCl with R = 0.4 at 5 Hz. 
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70 100 300 500 70 100 300 500 

A K (MPaVmm) A K (MPoN/mm) 

A) Forman Equation (closure). B) Forman Equation (no closure). 



A K (MPa'v/mm) 


C) SINK Equation. D) Sigmoidal Equation. 

Figure 22; Fit of fatigue data for Ti-6Al-4V in 1.0% NaCl with R = 0.7 at 5 Hz. 

66 


-I I 









A) Forman Equation (closure). 


B) Forman Equation (no closure). 


Figure 22a: Fit of fatigue data for Ti-6A1-4V in air with R = 0.1 at 5 Hz for p = q =2. 



A) Forman Equation (closure). 



100 300 500 700 

AK (MPaVmm) 


B) Forman Equation (no closure). 


Figure 22b; Fit of fatigue data for Ti-6A1-4V in air with R = 0.4 at 5 Hz for p = q =2. 
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70 100 300 500 

AK (MPaV mm) 


A) Forman Equation (closure). 


B) Forman Equation (no closure). 


Figure 22c: Fit of fatigue data for Ti-6A1-4V in air with R = 0.7 at 5 Hz for p = q = 2. 



A) Forman Equation (closure). 



B) Forman Equation (no closure). 


Figure 22d: Fit of fatigue data for Ti-6A1-4V in 1.0% NaCl with R = 0.1 at 5 Hz for p = 







A) Forman Equation (closure). B) Forman Equation (no closure). 


Figure 22e: Fit of fatigue data for Ti-6A1-4V in 1.0% NaCl with R = 0.4 at 5 Hz for p = 
q =2. 
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AK (MPaVmm) AK (MPaVmm) 


A) Forman Equation (closure). B) Forman Equation (no closure). 

Figure 22f: Fit of fatigue data for Ti-6A1-4V in 1.0% NaCl with R = 0.7 at 5 Hz for p = 
q =2. 
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Table II. 9: Coefficients of determination 
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Table II. 10 : Average coefficients of detennination. 


Equation Type 

Parameters Fit 

Environment 

Coefficient 

Forman (with closure) 

C, n 

Air 

0.8951 


(p = 0.25, q = 0.75) 

1.0%NaCl 

0.8816 

Forman (with closure) 

C, n, p, q 

Air 

0.9857 



1 . 0 % NaCl 

0.9777 

Forman (no closure) 

C, n 

Air 

0.8951 


(p = 0.25, q = 0.75) 

1 . 0 % NaCl 

0.8816 

Forman (no closure) 

C, n, p, q 

Air 

0.9857 



1 . 0 % NaCl 

0.9777 

Sigmoidal 

B, P. Q, D 

Air 

0.9857 



1 . 0 % NaCl 

0.9765 

- 

Hyperbolic Sine 

C 2 , C 3 , C 4 

Air 

0.9787 


{C, =0.5) 

1 . 0 % NaCl 

0.9459 


68 











































Table II. 10a: Coefficients of determination for p = q = 2.0 


Environment 

R 

Forman Closure 

Forman No Closure: 



C, n (p = q = 2.0) 

C, n (p = q = 2.0) 

Air 

0.1 

0.9034 

0.9034 

Air 

0.4 

0.8801 

0.8801 

Air 

0.7 

0.9846 

0.9846 

1.0% NaCl 

0.1 

0.9573 

0.9573 

1.0% NaCl 

0.4 

0.9128 

0.9128 

1.0% NaCl 

0.7 

0.7445 

0.7445 


Table II. 10b: Average coefficients of determination for p = q = 2.0 


Equation Type 

Parameters Fit 

Environment 

Coefficient 

Forman (with closure) 

C,n (p = q = 2.0) 

Air 

0.9227 

1.0%NaCI 

0.8715 

Forman (no closure) 

C.n (p = q = 2.0) 

Air 

0.9227 

1.0% NaCl 

0.8715 
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parameters were set equal to the appropriate linear relationship. In reality, only one 
data set is needed to determine the effect of frequency, however, data for 1 Hz and 10 
Hz were used so that the relationships determined would be accurate over a wider range 
of data than if only one data set was used. (If only one data set was used, the 
relationships determined would only be accurate for frequencies very near the frequency 
for which data were entered.) Since two data sets were used, the equation parameters 
were set equal to relationships involving two loading variables, frequency and either the 
hold time or the stress ratio, even though the hold time and stress ratio are constant. 
Using these relationships, the equation parameters were interpolated for a frequency of 5 
Hz. The fatigue crack growth was then calculated for a range of AK values using the 
resultant equation, and the interpolated growth rates were compared to actual data. 

The Forman Equation with closure was tested first, fltting only the parameters C 
and n; p was entered as 0.25 and q was 0.75. The fitted values are listed in Table II.ll. 


Table II.ll: Forman Equation with closure (C,n Fit) 


Parameter 

Environment 

1 Hz 

10 Hz 

c 

3.5% NaCl 

0.4398 X 10^ 

0.4007 X 10’ 

n 

3.5%NaCl 

3.481 

2.432 


The parameters were set equal to the following linear relationships: 


C=-0. 13741 JT10-® 



+ 0.45761A'10‘’' 


log 


X + 1 ' 

0.1 


(11.28) 
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(11.29) 


n = 4.040 



+ 1 


+ 2 . 265 log 


T + 1 

0.1 


= 13.648 log 


m 


+ 1 


+ 6.708 log 


T + 1 ' 

0.1 


(11.30) 


These results are graphically shown graphically in Figure 23. Next, the parameters C,n,p, 
and q were fit, and the resulting values are shown in Table 11.12. 


Table 11.12 Forman Equation with closure (C,n,p,qFit) 


Parameter 

Environment 

1 Hz 

10 Hz 

C 

3.5% NaCl 

0.5994X 10^ 

0.3430X 10^ 

n 

3.5% NaCl 

0.9971 

0.2544 

p 

3.5% NaCl 

2.971 

2.253 

q 

3.5% NaCl 

0.4662 

2.075 


The following relationships were determined: 


C=0.9877Z lO''* log 


B) 


+ 1 


+ 0.302122T 10 log 


T +1 ' 

0.1 


(11.31) 


n = 2 . 861 log 


(41 + l] + 0.1360 log 

A fJ J ^ [ 0 . 1 


(11.32) 


p = 2.765 log 1 + 2.138 log 


(11.33) 


q = -6.197 


log [(i) 


1 


+ 2.332 


log 

I 0.1 


(11.34) 
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da/dN (in/cycle) da/dN (in/cycle) 




itted data for 1 Hz. 


B) Fitted data for 10 Hz. 



5 7 10 30 50 70 100 


AK (ksi\/in) 



C) Interpolated data for 5 Hz. 


D) Effect of frequency on equation. 


Figure 23: Fitted and interpolated data using Forman Equation with closure (C,n fit) 







(11.35) 


Aj:,^ = 13.648 log + 1 


+ 6.708 log 


T + 1 ' 

0.1 


These results are shown graphically in Figure 24. 

The second equation tested was the Forman Equation without closure. Originally, 
only the parameters C and n were fit. The parameter m was held constant at 1.0, p was 
constant at 0.25 , and q was constant at 0.75. The fitted parameters are listed in Table 
n.l3. 


Table 11.13: Forman Equation without closure (C,n Fit) 


Parameter 

Environment 

1 Hz 

10 Hz 

c 

3.5%NaCI 

0.2008 X 10^ 

0.2335 X 10^ 

n 

3.5%NaCl 

3.481 

2.432 


The parameters were set equal to the following linear relationships: 


C=-0 . 82216 JT 10"'^ ^ + 0.267 57 ;sr lO’”^ log 


T + 1 ' 

0.1 


(11.36) 


72 = 4 . 040 log 


\ -f/ 


+ 1 


+ 2 . 265 log 


T + 1 


0.1 


(11.37) 


AJTfv = 13.648 log 


l—\ 
A W 


+ 1 


+ 6 .708 log 


T + 1 ' 

0.1 


(11.38) 


These results are graphically shown in Figure 25. Next, the parameters C,n,p, and q 
were fit, and the selected values are shown in Table II. 14. Again the parameter m was 
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A) Fitted data for 1 Hz. 


B) Fitted data for 10 Hz. 




O Interpolated data for 5 Hz. 


D) Effect of frequency on equation. 


Figure 24: Fitted and interpolated 


data using Forman Equation with closure (C,n,p,q fit) 
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held constant at 1.0. 


Table 11.14: Forman Equation without Closure (C,n,p,qFit) 


Parameter 

Environment 

1 Hz 

10 Hz 

c 

3.5%NaCl 

0.5028 X 10"* 

0.2887X 10^ 

n 

3.5%NaCl 

0.9971 

0.2544 

p 

3.5% NaCl 

2.971 

2.253 

q 

3.5% NaCl 

0.4662 

2.075 


The following relationships were determined: 


C=0 . 82467 JT 10'^ log 


U/ 


+ 1 


+ 0.25452.jr 0-^ log 


T + 1' 


0 . 1 


n = 2.861 log + ij + 0.1360 log 


p = 2.765 log + l] + 2.138 log 


T + 1 


0.1 


g = -6.197 log 1 + 2.332 log 


T + 1 


0.1 


= 13.648 log 


A W 


+ 1 


+ 6 .708 log 


T + 1 ' 

0.1 


(11.39) 


(11.40) 


(11.41) 


(11.42) 


(11.43) 


The results are depicted in Figure 26. 
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itted data for 1 Hz. 


B) Fitted data for 10 Hz. 




AK (ksiVin) 


C) Interpolated data for 5 Hz. 


Figure 26: Fitted and interpolated data using Forman Equation without closure (C,n,p,q 







The third equation tested was the Hyperbolic Sine Equation (SINK Equation). 
The values selected by the program are listed in Table n. 15, and results are graphically 
shown in Figure 27. 


Table 11.15: SINH Equation Parameters 


Parameter 

Environment 

1 Hz 

10 Hz 

Q 

3.5%NaCI 

4.9126 

4.3797 

C3 

3.5% NaCl 

-1.5793 

-1.4022 

C4 

3.5% NaCl 

-3.0435 

-4.0799 


The parameter Cj was held constant at 0.5. The following relationships were 
determined: 


Cj = 2 . 0527 log 


yj.\ 

1 fl 


+ l] + 4.3135 log \- 

' ' 0.1 


(11.44) 


C 3 = -0.6824 log 


(^) + l] - 1.3800 log \- O^l^l 

\fl J [0.1 


(11.45) 


Q = 3.9917 log 


(■i) * 1 ]- 4.2637 log[i-^ 


li? 


(11.46) 


The final equation tested was the Sigmoidal Equation. The values calculated by 
the program are listed in Table 11.16, and the results are graphically shown in Figure 28. 
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Tj-6A1-4V 
R = 0.1 

Frequency = 10 Hz 
3 . 5 % NoCI Solution 



AK (ksivin) 

C) Inteq)olated data for 5 Hz. 


AK (ksivin) 

D) Effect of frequency on equation. 


Figure 27: Fitted and interpolated data using Hyperbolic Sine (SINH) Equation. 







Table 11.16: Sigmoidal Equation Parameters 



The following relationships were determined: 

B = 2.3025 log j-ij ■ 8-0183 log ~ ~q ^ ^ (11.47) 

P = 4.0048 log [(41 + l] - 1.6660 log [ -^ ~ O.li? ] (11,48) 

\ff J [0.1 

Q = 2.6135 log + 1 + 2.1430 log ^ (11.49) 

D = 6.0157 log l] ~ 2.3375 log ^ (11.50) 

= 13 .64 8 log J-|| + 1 - 6.738 log ^ (11.51) 

In all of the tests, the value of Kc was 83.27 ksi/in. The values of AK* were entered as 
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A) Fitted data for 1 Hz. 


B) Fitted data for 10 Hz. 




Q Interpolated data for 5 Hz. D) Effect of frequency on equation. 

Figure 28: Fitted and interpolated data using Sigmoidal Equation. 
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10.82 ksiv/in for 1 Hz and 7.27 ksi/in for 10 Hz. The value of AK«, was interpolated for 
5 Hz. 

As seen in Figures 23-28, the actual data for T1-6A1-4V in 3.5% NaCl solution 
display the frequency "crossover" effect first reported in a -I- |8 titanium alloys by Dawson 
and Pelloux“‘, and later by Dawson^® and Yoder and coworkers^^. The "crossover" effect 
is a change in the ordering of da/dN values as a function of frequency at a transition AK 
value associated with the onset of "cyclic" stress corrosion cracking (SCC)*“. Dawson 
and Pelloux called this transition level AKjcc. a frequency dependent parameter lower 
than Kjscc- Below the transition AK level, da/dN decreases with decreasing frequency 
due to passive film formation at the crack tip‘‘‘. Above AKjcc. da/dN is controlled by 
hydrogen embrittlement mechanisms, thus da/dN increases as frequency decreases and 
more time is allowed for hydrogen embrittlement to occur^’. The interpolated data also 
displayed this "crossover" effect. 

Tests were also conducted in which the data for 1 Hz and 5 Hz were entered, and 
data for 10 Hz were interpolated, as well as tests where the data for 5 Hz and 10 Hz 
were entered, and data for 1 Hz were interpolated. Thus, three cases were analyzed: 

1) Case 1: the frequency at which data were interpolated was lower than the 

frequencies for which data were entered 

2) Case 2: the frequency at which data were interpolated was between the 

frequencies at which data were entered 

3) Case 3: the frequency at which data were interpolated was higher than the 

frequencies at which data were entered 
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The best results occurred for Case 2, where the frequency at which data were 
interpolated is between the range of frequencies for which data were entered. This is 
the case for which the results were presented in this report. When the frequency at 
which data were interpolated was higher than the frequencies at which data were 
entered, the results were less accurate. When the frequency at which data were 
interpolated was lower than the frequencies at which data were entered, the results were 
least accurate based upon the coefficients of determination. 

It should be noted that the ability of the model to interpolate fatigue behavior in 
material-environment systems where the shape of the da/dN versus AK curve is strongly 
dependent on the environment has not been determined. Instead, the curves studied 
followed the commonly analyzed air case in shape. 


n.4: Conclusions for the Interpolative Model 

One aspect of this project was to develop an interpolative model for use by NASA 
in analyzing and predicting the fatigue behavior of metals in corrosive environments. A 
stand-alone computer program based on NASA FLAGRO was developed which can be 
used to curve fit da/dN versus AK data as well as interpolate trends in fatigue behavior. 

1. The regression subroutines used in the linear least squares analysis of 

the Paris, Forman, and Sigmoidal Equations were bench marked against the 
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curve fitting subroutines contained in the FLAGRO program. Both the 
stand alone program and the FLAGRO program generated the same 
results for C and n when curve fitting a data set to the Forman Equation 
with closure. 

The non-linear least squares algorithm was bench marked against results 
published by Haritos*®. Using the baseline data from Haritos’® , the values 
for the parameters in the SINH Equation were optimized. The values for 
the parameters determined by the interpolative program were different 
than the results published by Haritos^®, however, the results generated by 
the stand-alone program yielded a slightly better fit based upon the 
coefficients of determination. 

Based on the data for Ti-6A1-4V in air and a 1 % NaCl solution, the 
Forman Equation (with and without closure) yielded the best results when 
fitting C, n, p, and q, with an average coefficient of determination of 0.9817 
for all data tested. The Sigmoidal Equation displayed similar accuracy, 
with an average coefficient of determination of 0.9811 for all data tested. 
The SINH Equation also fit the data with accuracy, yielding an average 
coefficient of determination of 0.9623. The Forman Equation (with and 
without closure) was least accurate when the parameters C and n were fit; 
the coefficient of determination was 0.8884. 

All of the crack growth rate equations were able to interpolate trends in 
da/dN versus AK behavior when the frequency at which data were 



inteipolated was within the range for which data were entered. The model 
was not able to accurately predict data when the frequency at which data 
were interpolated was outside the range at which data were entered. 

The data for Ti-6AI-4V in 3.5% NaCl solution, used to test the model, 
displayed the frequency "crossover" effect. The interpolated data also 
displayed this effect. The model accurately shows that AK* increases as 
frequency decreases. 

Future work will concentrate on improving the logarithmic assumption 
which relates the equation parameters to the loading characteristics. 

This will be done by generating corrosion fatigue data for Ti-6A1-4V (ELI) 
in an aqueous chloride environment. These data will then be used to 
improve the relationships by analyzing the effects of frequency, hold time, 
and stress ratio on (da/dN)tou>j- 



Section III: Linear Superposition Model 


in.l: Introduction to the Linear Superposition Method 


Although there are several methods of predicting environment-enhanced fatigue 
crack growth in materials, the simplest is linear superposition. This method predicts 
crack growth by summing the effects of mechanical fatigue and stress corrosion. Wei 
and Landes first proposed that the effect of an aggressive environment on fatigue crack 
propagation in high strength steels is independent of cyclic loading. Both effects can be 
added to give a first-order approximation of the total environment enhanced crack 
growth rate^. In equation form; 


da _ da ^ da 

total fatigue stzeBa corrosion 


(III.l) 


Since stress corrosion crack growth is measured per unit time, it must be calculated for 
an entire cycle of fatigue loading. If the stress corrosion crack velocity, (da/dt)j„yjronn,en[, 
can be expressed as a function of stress intensity, K, this equation can be integrated over 
a cycle of loading to create (da/dN),^^, In equation form: 


da 

dN 

8tie39 coxTosion 



[jr( t) ] dt 


(III. 2) 


where t is the period of the fatigue load cycle. In actuality, environment may only affect 
cracking during a portion of this time; for example during loading but not unloading. As 
an example, if the crack growth velocity is a constant rate for the range of stress 
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intensities in a fatigue load cycle: 


^ = C for all K(t) (III. 3) 

at 


then crack growth per cycle is given as: 

BtiesB coixosion 



C T 


(III. 4) 


and depends only on the period of the loading cycle, t (t = 1/f), and the constant C. 
For a more complex relationship between crack growth velocity and K such as: 

= CK{t)° (III. 5) 

dt 


the stress corrosion crack growth per cycle depends on the period and the loading 
function as well as the constants C and n. As shown in Appendix E, solving Equation 
in.2 in conjimction with Equation HI. 5 with n = 2 and sinusoidally varying stress 
intensity yields: 

da 

tress cozzosBlon 

where K^vg = + Kn„n) ! 2. Additional solutions are derived in Appendix E for 

square wave and ramp loading coupled with a variety of n values. 

The success of this method of prediction has been conErmed for various material- 
environment combinations^-^*'^’. When Wei aiKl Landes originally proposed the linear 
superposition method, they supported their work with fatigue crack growth rate 
predictions for 18Ni (250) maraging steel in gaseous hydrogen. The fatigue component 


c X 




+ K‘ 


avffl 


(III. 6) 


of the prediction used data for an 
argon environment which was 
frequency independent and 
representative of purely mechanical 
fatigue damage, as seen in Figure 29. 

Constant load time-based crack growth 
rate data for 18Ni steel in 
dehumidified hydrogen were 
integrated over a cycle of loading to 
create the stress corrosion component 
added to the baseline fatigue data^. 

The procedure reasonably predicted the absolute magnitude and frequency dependence 
of H 2 -enhanced fatigue crack propagation in 18Ni (250) maraging steel as shown in 
Figure 29. i 

The effects of AK and frequency on crack growth rates for AA 7079-T651 in 
saturated NaCl solution were predicted by Speidel using thenmethod of superposition^®. 
His predictions, displayed in Figures 30 and 31, show the accuracy of the superposition 
model, especially at frequencies less than 1 Hz. 

Harmon and coworkers” employed the linear superposition method to predict the 
effects of AK, hold time, and loading frequency on elevated temperature: fatigue crack 
growth rates in nickel-base superalloys. This approach was added to an existing fatigue 
crack growth life prediction program, CRKGRO. To implement this approach, the Wei- 



te/iN. H. FCR CVCLE 


Figure 29: Measured and linear 
superposition predicted fatigue 
crack growth rates for maraging 
steel in H,. After Wei and Landes^. 
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Landes model was used to predict fatigue 
crack growth in IN718 and AF115 at 
1200” F. For both examples, shown in 
Figures 32 and 33, the effects of hold time 
are shown to be effectively modeled by 
superposition^’. 

The preceding examples show the 
capability of the linear superposition 
method to predict fatigue crack growth 
rates for structural alloys in various gas 
and liquid environments, including 
elevated temperature. The superposition 
method is, however, only accurate and useful for limited cases. Specifically, the alloy 
must be extremely sensitive to stress corrosion cracking. The contribution of stress 
corrosion to crack growth must be significantly greater than that of inert environment 
mechanical fatigue. This condition is achieved when the loading frequency is low, the 
stress intensity is much higher than Kiscc» ^nd da/dt is relatively rapid. The predictions 
made by Speidel were for AA 7079-T65 1 , a material which is affected strongly by 
environment compared to other 7xxx aluminum alloys, making it a good candidate for 
this modeP. The maraging steel is similarly sensitive to gaseous hydrogen embrittlement 
and nickel-based superalloys are prone to sustained load crack growth at relatively low 
applied K. Linear superposition is not likely to be useful for most structural alloy- 



linear superposition predicted 
effect of loading frequency on 
NaCl enhanced da/dN„,^] for 
AA7079. After Speidel^*. 
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Figure 31: Linear superposition predicted (solid lines) and 
measured (data points) fatigue crack growth for stress corrosion 
cracking prone AA7079 in NaCl . After Speidel“. 
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Figure 32 Application of Wei-Landes 
linear superposition model to fatigue 
crack propagation in IN718 at 1200" F 
with varying AK,R and t. After 
Harmon^’ . 



Figure 33: Application of linear 
superposition to predict the effect 
of hold time (t) on da/dN for IN718 
at 1200° F. After Harmon^’. 
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environment systems where the material is resistant to stress corrosion cracking. 


m.2: Implementation of the Linear Superposition Model 

The linear superposition method has been implemented in a FORTRAN program 
to predict corrosion fatigue crack propagation rates. The program offers the user 
options in determining the variables which affect crack growth rates. These options 
include; 1) the equations used to represent inert environment mechanical fatigue and 
stress corrosion crack growth rates, 2) the type of loading to which the cracked 
component is exposed, 3) the period of the loading cycle and the portion of this period 
over which environmental cracking occurs, 4) the range of output stress intensity factor 
ranges (AK), and 5) the number of output data points. With these inputs, the program 
predicts the total fatigue crack growth per cycle of loading for a given material- 
environment system. The method predicts the effects of AK, stress ratio (R), frequency 
(f), hold time (t) and loading wave form. 

l ^epresentiriQ Crack Growth in Equation Form 

The linear superposition method assumes that crack growth due to stress 
corrosion and inert environment mechanical fatigue are independent of one another. 

Both of these independent quantities are calculated over one cycle of loading and are 
added to give (da/dN)toai- The program allows for four different methods of 
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representing the crack growth rate laws. 


A) Single Power Law : The single power law, or Paris Equation, is the most 
common method of representing mechanically-based fatigue crack growth 
rates. This method assumes that a linear relationship exists between the 
logarithms of the growth rate and the stress intensity factor range. In 
equation form: 


da 

dN 


AAiif® 


(III. 7) 


An equation of this form can also be used for stress corrosion crack growth 
by interchanging da/dt for da/dN and replacing AK with K. 

B) Double Power Law: The double power law uses the same relation between 
stress intensity and crack growth as the single power law (Equation HI. 7) 
except the constants A and m change when the stress intensity reaches a 
certain level. By using two Paris Equations, the program can more 
accurately represent fatigue data near threshold. A two power-law 
approach is also useful for representing da/dt versus applied K. 

C) Forman Equation 

D) Sigmoidal Equation: To utilize either the Forman or the Sigmoidal 
Equations, which were discussed in the interpolative modeling section, the 
program must establish the equation constants for a specific material. The 
program allows the user the option of entering in the values of these 
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constants, or inputing a set of data and allowing the program to calculate 
the constants by linear regression analysis. When using the Forman 
Equation, the user is offered the option to read the constants from a 
library of materials that was extracted from the NASA-FLAGRO program. 

For the program to function properly, the input da/dN and AK data must be in 
two columns. Whenever data are input from a file, the program assumes the first 
column contains the stress intensity factor range values, AK. The program allows the 
user to alter data files in three ways. An option is available to add an end of file (EOF) 
marker to the data. Most data files contain EOF markers when created, however this 
option insures that the program will not return an error while importing data from a file. 
The sorting option allows the user to order the data from smallest to largest stress 
intensity values. The conversion option allows the user to multiply each column of data 
by a constant factor as a method of converting units. 

For any situation in which data are used to determine equation constants, the 
program allows the user to view a graph of the data along with the fitted equation. Each 
axis can either be spaced linearly or logarithmically. Points from the equation, generated 
in the plotting of the curve, can also be written to a file. 

Representing the Load Wave Form 

A wide variety of time dependent loading forms are important to components in 
service. For simplicity, only a few loading forms are typically employed in laboratory 
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characterization of corrosion fatigue. This program allows the user to choose from three 
different types of loading. Additional forms can be added if necessary. 

A) Sinusoidal Loading: This method of loading is a standard experimental 
practice, with the load constantly varying with no lu)ld time. The stress 
intensity is modeled by: 


Kit) = 


2 


COS 


2ti t 

X 



(III. 8) 


where t is the period of the load cycle. The load cycle begins at K„i„, the 
applied stress intensity increases to a maximum (K,^ at t/2, then 
decreases to Kn,j„ at the end of the period as seen in Figure 34. 

B) Ramp Loading: Ramp loading allows the user to control the percentages 
of the load cycle involved with increasing K from K„,i„ to and with 
maintaining a constant K-hold period at By default, the program 
assigns the remaining percentage to the unloading portion of the fatigue 
cycle, as seen in Figure 35. This function allows for no hold time, if 
requested, to simulate sawtooth or triangular loadings, as seen in Figure 36. 

C) Square Wave Loading: Square wave loading allows for periods of constant 
maximum and minimum stress intensities with negligible loading and 
unloading times. The program allows the user to select the percentage of 
the time spent at maximum stress intensity and allots the remaining time in 
the period to constant minimum stress intensity. A typical example is 
shown in Figure 37. 
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Stress intensity 



constant K^. 
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Figure 36; Ramp loading with no hold time. 


Square Load Waveform 



time 

Figure 37; A typical square wave loading function. 
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Percent Reduction in Loadins 


While results are not conclusive and general for all material-environment systems, 
it is likely that stress corrosion cracking only contributes to the environmental da/dN 
during the loading and constant K hold-time portions of a load cycle^. This notion is 
reasonable because initial unloading of a sharp crack causes a compressive stress 
singularity local to the tip. Environmental cracking is only operative under tensile 
stresses which are produced when the cycle again reverses to the loading mode. The 
program, therefore, allows the user to input a percent reduction on unloading crack 
growth. Initially, the program calculates the crack growth contribution from stress 
corrosion during unloading equivalent to that for loading. The input reduction is then 
taken into account by eliminating a percentage of the stress corrosion crack growth rate 
calculated for unloading. 


Output Data Parameters 

The program creates an output file which contains predicted corrosion fatigue, 
da/dNy,^, versus AK values as well as additional information. The user has the following 
options in controlling the output data. 

A) Constant R, K„^, or K„,jn; The output data file contains two columns of 
data; the X column of stress intensity ranges, AK, and the corresponding 
Y column of corrosion fatigue crack growth rates, (da/dN)^. Since 
fatigue crack growth rates are a function of the stress ratio, R, a ratio must 
be specified for each data point for the file to be useful. Most commonly, the 
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stress ratio is constant for a data set; however, if a constant or is 

specified, the stress ratio can be calculated for each AK and the input 
constant or 

B) Maximum and minimum AK values. 

C) Number of data pairs. 

D) Spacing of AK values: This option allows the user, if predicting (da/dN)t(,aj 
over a large range of AK values, to space the AK data evenly on a linear or 
logarithmic scale. 

D) Frequency of loading: For all loading conditions, the program will ask for 
the frequency, f, of the loading cycle. The reciprocal of the frequency is 
the period, r, therefore, the user must either know the loading frequency 
of interest, or know the time required to complete a cycle and take the 
reciprocal as the frequency. 

The Calculation - 

Once the user has input the required parameters, the program calculates 
corrosion fatigue crack growth rates. For each AK specified by the user to be output, 
the program first calculates the crack growth rate due to fatigue. Hiis is calculated using 
the equation chosen by the user to model fatigue. Next, the stress coSbsion'compbneht 
is calculated by integrating the equation used to model time-based stress corrosion. For 
each loading condition, the integration is different. With sinusoidal loading, crack 
growth is integrated over the loading portion of the cycle (period = l/2f). The 
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unloading crack growth is initially equated to the loading portion, then adjusted for the 
input percent reduction. With ramp loading, crack growth during increased load is 
calculated first, followed by a calculation of crack growth due to the static load during 
the hold time. The unloading contribution is integrated, then adjusted for the percent 
reduction. For square loading, the calculation sums the stress corrosion for two different 
K values and K„,in) over a period of loading. No reduction is calculated since the 
time period of unloading is zero. Once the crack growth rate due to stress corrosion for 
each portion of loading is calculated, the fatigue component is added to predict the 
corrosion fatigue crack propagation rate. 


m.3: Data Testing and Results for the Linear Superposition Model 

As shown earlier, the linear superposition method was used by Wei and Landes^, 
SpeideP, and Harmon^^ to successfully predict fatigue crack growth rates in corrosive 
environments. To prove the capability for reliable predictions, the current program was 
first tested to insure that the integration method was accurate. The program was then 
employed for the same material-environment combinations tested by Wei and Landes^, 
SpeideP, and Harmon^’ in an attempt to duplicate their results. 

The effort to test the accuracy of the program focused on the numerical 
integration component. Since the superposition model works best when predicting crack 
growth rates with a majority contribution due to stress corrosion, inaccuracy in the 
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integration of stress corrosion over a loading cycle would cause significant errors in the 
final results. To test the integration, the exact solutions of hypothetical integrations were 
compared to program calculations. The exact solutions for all loading conditions were 
calculated for various models of stress corrosion crack growth. As shown in Appendix E, 
program results were completely accurate aside from round-off errors. 

To obtain the baseline fatigue and stress corrosion models used in the original 
predictions, either the equation constants for one of the models or the original 
experimental data were required. In most cases only plots of data and modeled 
predictions were available from the original reports. To obtain the necessary 
information, graphs of data taken from the reports were digitized. 

The computer program was first compared to the linear superposition calculations 
in the Wei and Landes paper^. The base fatigue data used were for 18Ni (250) maraging 
steel tested in argon. This was modeled by a two-slope Paris Equation. The stress 
corrosion data, also obtained from the Wei and Landes paper^, were modeled by a 
Sigmoidal Equation. As Figures 38 and 39 show, the program’s models of the corrosion 
fatigue crack growth rates are consistent with the actual data, although the predictions do 
not coincide with the original predictions. Two differences in the calculations may 
explain this discrepancy. First, the program used Paris Equations to model the fatigue 
data and a Sigmoidal Equation to model stress corrosion kinetics. The original 
predictions made no transformation of data to an equation. Second, the original 
predictions were made by approximating the integral of stress corrosion over one cycle of 
loading^. Since Wei and Landes only meant for the superposition model to be a first- 
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Figure 38: Linear superposition predictions, compared to measured (da/dN)uKai. for 
18Ni(250) maraging steel in gaseous hydrogen. 
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Figure 39; Linear superposition predictions, compared to measured (da/dN)uRai. for 
18Ni(250) maraging steel in at a lower loading frequency. 
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order approximation of crack growth rates, the need for precision was minimal. The 
current program uses Simpson’s Rule to more accurately evaluate the integral. 

The computer program results were then compared to predictions by Speidel, 
using the linear superposition model for AA7079-T651 in saturated NaCl solution, and 
focusing on the strong effect of loading frequency^. To simulate Speidel’s predictions, 
the stress corrosion data were approximated by a double power law. Speidel used a 
fatigue equation including the critical stress intensity value for unstable fracture, Kc, and 
the threshold stress intensity range, AK^, to model (da/dN)faUgue • The program does not 
use this equation to model fatigue; however, the Forman Equation also incorporates the 
critical and threshold values for a material. Since the program will fit the Forman 
Equation to an input set of data, several points along the fatigue model Speidel used 
were calculated and fit to the Forman Equation. This approximation was sufficient since 
Speidel ignored the negligible fatigue contribution to crack growth at low frequencies^®. 

The results from the program are comparable to those reported by Speidel as 
shown in Figures 40 and 41. The difference between the prediction and measured 
results at 1 to 10 Hz is attributed to the effect of time-cycle dependent corrosion fatigue 
which is not accounted for in the superposition model. 

A third comparison was made to predictions published by Ashbaugh and 
Nicholas^*. They observed the temperature dependence of sustained load cracking in 
IN718 at elevated temperature and proposed that the superposition model would work 
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Figure 40: Hie loading frequency effect on the total environment-enhanced fatigue crack 
growth rate for AA7079 in aqueous chloride. Measurements reported by SpeideP are 
compared to the solid line predicted by the linear superposition model. 
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Figure 41: Linear superposition predictions of corrosion fatigue crack growth 
kinetics in a stress corrosion cracking susceptible aluminum alloy. 
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well in predicting the effect of temperature on fatigue crack growth rates. Ashbaugh and 
Nicholas^’ modeled both the fatigue crack growth and stress corrosion with a Sigmoidal 
Equation. The predictions for FCP at two different temperatures were promising. 

The program’s predictions correspond well with those of Ashbaugh and Nicholas^’, 
as seen in Figure 42. The loading frequency and wave form were duplicated from the 
experimental description. Sustained load cracking was not considered to occur during 
the unloading portion of the load cycle (crack growth period = l/[2f]). For stress ratios 
of 0.1 and 0.5 the predicted growth rates agreed with measured values. 


in.4: Conclusions for the Linear Superposition Model 

A computer program was created to predict corrosion fatigue crack propagation 
rates for materials in corrosive environments using the linear superposition method. 

1. Although linear superposition is not accurate for all material-enviro nm ent 
systems, this method has proven to be reliable in those limited cases 
where the environment has a strong effect on subcritical crack growth. 

2. The program incorporates conditions which affect environmental fatigue 
cracking, including AK, R, f, t, wave form and material fatigue and stress 
corrosion behavior. 

3. The integration portion of the program was tested and proven to accurately 
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Figure 42: Linear superposition predictions of the effect of stress ratio on elevated 
temperature environmental fatigue cracking in IN718. Measured data were reported 
by Ashbaugh and Nicholas^*. 




calculate the total crack growth due to stress corrosion over a cycle of 
loading. 

4. Comparisons to various applications of superposition verify the accuracy 
and capability of the current program to model experimental fatigue data. 
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Appendix A: Data File Directory 


4340 Steel 


File 

X- 

Frequency fHzl 

Environment 

Units 

T41A 

0 

0.1 

Air 

m/cyc, Pa-s/m 

T41B 

0 

0.1 

3.5% NaCl 

m/cyc, Pa>/m 

T41C 

0 

0.1 

Dist. water 

nu/cyc, Pa^^l 

T42A 

0 

0.001 

Dist. water 

m/cyc, MPaVm 

T42B 

0 

0.01 

Dist. water 

m/cyc, MPa>^m 

T42C 

0 

0.1 

Dist. water 

m/cyc, MPa-v^m 

T42D 

0 

1 

Dist. water 

m/cyc, MPaVm 

T42E 

0 

4 

Dist. water 

m/cyc, MPa-vAm 

T42F 

0 

4 

Vacuum 

m/cyc, MPa-v/hi 

T43A 

0 

0.4 

Argon 

m/cyc, MPa>/m 

T43B 

0 

0.4 

Dist. water 

m/cyc, MPav^i 

T43C 

0 

0.2 

Dist. water 

m/cyc, MPa-/m 

T43D 

0 

0.04 

Dist. water 

m/cyc, MPa-sAm 

T44A 



3.5% NaCl 

in/min, ksi-s/in 
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Ti - 6A1 - 4V 


File 

R 

Frequency 

TTIA 

0.1 

10 

TTIB 

0.1 

10 

me 

0.1 

5 

TTID 

0.1 

1 

TT2A 

0.5 

2 

TT2B 

0.5 

2 

TT2C 

0.5 

2 

TT2D 



TT2E 



TT3A 

0 

4 

TT3B 

0 

4 

TT3C 

0 

4 

TT4A 



TT5A 

0.05 

50 

TT5B 

0.05 

50 

TT5C 

0.05 

50 

TT5D 

0.05 

50 


Environment 

Units 

Air 

in/cyc, ksi>/in 

3.5% NaCl 

in/cyc, ksiVin 

3.5% NaCl 

in/cyc, ksi-s/in 

3.5% NaCl 

in/cyc, ksiVin 

3.5% NaCl 

cm/cyc, ksiVin 

Dist. water 

cm/cyc, ksiVin 

Argon 

cm/cyc, ksi-v/in 

Dist. water 

cm/sec, ksiVin 

3.5% NaCl 

cm/sec, ksiVin 

Argon 

cm/cyc, ksiVin 

Dist. water 

cm/cyc, ksiVin 

3.5% NaCl 

cm/cyc, ksi-v/in 

3.5% NaCl 

cm/sec, ksiVin 

Argon 

in/cyc, ksiVin 

Air 

in/cyc, ksiVin 

Dist. water 

in/cyc, ksi-y/in 

3.5% NaCl 

in/cyc, ksiVin 
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7075 Aluminum 


File 

T71A 

T71B 

T72A 

T72B 

T72C 

T72D 

T72E 

T72F 

T73A 

T73B 

T73C 

T73D 

T73E 

T73F 

T73G 



Frequency (Hz) 

Environment 

Units 

0.1 

5 

Air 

MPa-^m, m/cycle 

0.1 

5 

Argon 

MPa-v/m, m/cycle 

0 


Air 

MPa-^/m, mm/cycle 

0.2 


Air 

MPa-v/m, mm/cycle 

0.33 


Air 

MPa-N^m, mm/cycle 

0.5 


Air 

MPa-v/m, mm/cycle 

0.7 


Air 

MPa-/m, mm/cycle 

0.8 


Air 

MPa-v/m, mm/cycle 

0.1 

5 

Dist. water 

MPa-v/m, m/cycle 

0.1 

5 

Air 

MPa-v/m, m/cycle 

0.1 

5 

H 2 O vapor (67 Pa) 

MPa-v/m, m/cycle 

0.1 

5 

H 2 O vapor (13 Pa) 

MPa-v/m, m/cycle 

0.1 

5 

HjO vapor (4.7 Pa) MPa-v/m, m/cycle 

0.1 

5 

HjO vapor (2.7 Pa) MPa-/m, m/cycle 

0.1 

5 

HjO vapor (2.0 Pa) MPa>/m, m/cycle 
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Appendix B: Procedure for Creating Data Files from Graphs 


I. DIGITABLET USE 

A. Obtain a UNIX account from the ACC. 

B. In the basement of the Albert H. Small building, turn on the Numonics 

Digitablet. The control is located on the back side of the console used to 
display the digits. 

C. Reset the digitablet on the back right comer of the table top. 

88888888 88888888 will appear on the console. 

D. Put the cross hair cursor over the cross hair at the bottom left comer of the 
menu and press MENU ORIGIN. 

E. Numerical values will now be shown on the console. Move the cursor and 

center the graph so that its axes are parallel to those on the digitablet (i.e. 

when the cursor is moved along the x-axis,a constant y value is shown on the 
console). 

F. Fasten the graph to the table top. 

G. If there are x and y values at the graph origin, place the cursor over the origin 
and press ’O’. If there are no labeled x and y values at the origin, put the 
cross hair over the intersection of the first x and y values on each axis. 

H. Select ENGLISH or METRIC units on the menu by pressing the READ 
MENU key on the cursor. 
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Record the characteristics of the graph on a reference sheet. The 
environment and any factors specific to the graph such as ftequency, pressure, 
or stress ratio should be noted. 

1. Record the units on each axis. 

2. Record the x and y values at the origin (as entered by the cursor). 

3. Another reference point on each axis is needed. The highest labeled 
value on each axis should serve as the reference point. Improved 
accuracy is achieved by choosing the point as far away from the origin 
as possible. Record the graph value of each point. 

4. For each axis, the distance between the origin and the reference point 
should be recorded. The distance can be read from the digitablet 
console. 

Use the cursor and the READ MENU key to select DIGITIZING. 

Logon to a network computer (not the computer Mxt to the digitablet). At 
the login prompt, type 

"telnet" <cr> (carriage return) 

Alt-A 

"kelvin.seas" <cr> 

Enter your UNIX login ID and then your password. Type <cr> when 
prompted for the terminal type. Keep in mind that capital letters have 
different meanings than lower case letters. 



Turn on the computer next to the Digitablet and type 
<cr> 

"connect telnet" <cr> 

"connect kelvin.seas" <cr> 

Enter the login ID and password again, 
type <cr> 

"export TERM " < cr > 

Now invoke the Jove editor along with the name of the data file that you are 
about to create. Type 

" j ove Filename . DAT " < cr > 

On the digitablet, enter points by placing the cross hair over the point and 
pressing ENTER OUTPUT. 

When all points have been digitized, enter the last point with negative x and 
y coordinates to signal the end of the data file for the transformation 
program. If negative values of points need to be digitized, the final point 
should be less than the lowest value. The process will function without the 
negative number at the end of the file, but it prevents any points from being 
inadvertently entered before the file is closed. 

Save the file and exit the editor. 

Control-X S to save 

Control-X Control-C to exit Jove 

If another set of data points off the same graph is desired, return to step P. 



If points from another graph are desired, repeat steps E through I before 
moving ahead to step P. 

Exit UNIX by typing "exit" <cr>on both terminals. If you wish to view the 
files that you have saved while in UNIX, type "Is" <cr>. 



II. REMOVING DATA FILES FROM UNIX 

A. The data files must be removed from UNIX and sent to a disk. While in 
UNIX, type "Is" <cr>to get a list of the files. 

B. To get into the file transfer mode, type 

Alt-T (initiates file transfer) 

Alt-I (send internet address) 

Alt-W (ask for password) 

Alt-W (send password) 

C. To transfer a file, type 

"put" <cr> 

At the local file prompt, type the UNIX "Filename. DAT" <cr>. At the 
distant filename prompt, use the drive you would like to send the file to 
{"B-.Filename.DAT"). 

D. To transfer another file, return to step C. 

E. Type "quit" <cr>to return to the UNIX operating mode. 

F. Type "exit" to leave UNIX. 
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in. DATA FILE TRANSFORMATION 


A. The program is available in an executable form on the program disk. To run 
the programs, type MENU. A title screen will appear, followed by a main 
menu. Select the option to digitize a set of data at the prompt. 

B. Follow the instructions and enter the data that were recorded during the 
digitization process. The data file that was retrieved from UNIX is not erased 
by this program. If a mistake is made, run the program again. 

C. If the program needs to be altered due to negative values, or any other 
reason, follow the steps below. 

1. At the F. /og/n prompt type "watfor77" <cr> 

2. Type "edit DIGITIZE.FOR" <cr> to retrieve the chosen FORTRAN 
program 

3. Use the function key F4 to go into and out of the program. 

4. Type "run" <cr> to start the program. 

5. "quit" <cr> exits FORTRAN without saving changes in dbie program, 
"exit" <cr> exits and saves the changes in the program. 


Appendix C: DIGITIZE Program Listing 

*$noex 

program DIGITIZE 
*$nolist 

*$include Gksdefn.for 
*$list 

* Jonathan M. Pope, Materials Science Department, University of 

* Virginia, June 1993. 

* 

* This program is to be used in conjunction with the Numonics 

* 2400 Digitablet. The Digitablet output can be transformed 

* into digital data representative of an original graph. Any 

* graph that has been digitized can have the values of its 

* points calculated by using this program and the data 

* requested in the digitization manual. Any combination of 

* linear and logarithmic scales on the X and Y axes can be 

* handled. 

* 

Real X1(1000),X2(1000),Y1(1000),Y2(1000) 

Real ORX,ORY,ONEX,ONEY,VX,VY 
Real LOGR,LINEAR 
Integer I,J,K,AX,AY 
Character*32 INPUT,OUTPUT 

* 

* INPUT is the file that was retrieved through the editor when 

* using the Digitablet. OUTPUT should be named by the user. 

* 

Call Screen 

print *, ’DIGITIZE Program - Jonathon Pope’ 
print *, ’University of Virginia - June 1993’ 
print *, ’Department of Mechanical Engineering’ 
print * 

print*, ’Enter the drive and name of the source data file’ 
read (*,100) INPUT 
1(X) format (a32) 
print* 

print*, ’Enter the drive and name of the desired output file’ 

read (*,100) OUTPUT 

print* 

print*, ’If the X axis has a logarithmic scale enter 1,’ 
print*, ’or enter 0 for a linear scale’ 
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read*, AX 
print* 

print*, ’Enter the X axis value at the origin’ 

read*, ORX 

print* 

print*, ’Enter the X axis value at the X axis reference point’ 

read*, VX 

print* 

print*, ’Enter the Digitablet console value for the distance’ 
print*, ’between the origin and the reference point’ 
read*, ONEX 

AX, AY are a logical operators used to remember if the axis 
is logarithmic or linear. ORX, ORY are the origin values on 
each axis. This is the origin as defined by the digitablet, 
and it may be different than the graph origin. VX, VY are 
each an additional referent point on each axis, arel ONEX, 
ONEY are the distances from the origin to the reference point. 
These distances are in centimeters, simply to be consistent 
with the manual, and are easily read from the console. 


print* 

print* 

print* 

print*, ’If the Y axis has a log scale enter 1,’ 
print*, ’or enter 0 for a linear scale’ 
read*, AY 
print* 

print*, ’Enter the Y axis value at the origin’ 

read*, ORY 

print* 

print*, ’Enter the Y axis value at the Y axis reference point’ 

read*, VY 

print* 

print*, ’Enter the Digitablet console value for the distance’ 

print*, ’between the origin and the reference point’ 

read*, ONEY 

open (l,file=INPUT) 

open (2,file= OUTPUT) 

1=1 

Xl(l)=.001 


The manual calls for the last point to be entered as a 
negative value off the graph. It is not necessary for the 



* program to function properly, but it is helpful if the file 

* needs to be investigated for errors. In the event that 

* negative valued points are part of the data file, the 

* parameter in the ’while’ statement below should be changed 

* to a smaller number. 

* 

while (Xia).ge.-1.0)do 
1 = 1+1 

read (l,*,end=400) Xl®, Y1(I) 
end while 
400 close (1) 

* 

* The mathematical functions needed to correct the data 

* for each axis are activated in the statements below. 

* 

J=I-1 

do 200 I=2,J 

if (AX.eq.l) then 

X2(I) = LOGR (Xia),ORX,ONEX,VX) 
else 

X2(I) = LINEAR (Xl(I),ORX,ONEX,VX) 
end if 

if (AY.eq.l) then 

Y2(I) = LOGR (Yl(I),ORY,ONEY,VY) 
else 

Y2(I) = LINEAR (Yl(I),ORY,ONEY,VY) 
end if 
200 continue 

* 

* The transformed data is written into the output file. 

* 

N = J - 1 

Call Sort (X2,Y2,N) 
do 300 K=2,J 

write (2,310) X2(K),Y2(K) 

310 format(G12.4,3X,G12.4) 

300 continue 

stop 

end 

* 

* This function corrects the digitablet values for axes 

* with logarithmic scales. 
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♦ 

real function LOGR(Z,OR,ONE,V) 
real Z,OR,ONE,V 

LOGR = OR*((V/OR)**(Z/ONE)) 

return 

end 


* This function corrects the digitablet values for axes 

* with linear scales. 

i|i ♦♦♦♦♦************i|ci|.******i|c*******Ktiii*i|.****j(.*****4,**,|,+**** 

* 

real function LINEAR(Z,OR,ONE,V) 
real Z,OR,ONE,V 

UNEAR = OR + (V-OR) * Z/ONE 

return 

end 

Subroutine Sort (X,Y,N) 

* 

* This subroutine sorts the Delta K and da/dN arrays so that the Delta 

* K values are placed in ascending order. 

Real X(N),Y(N), First, Tempx.Tempy 
Integer I, Index,! 

Do 30 I = 2,(N-1) 

First = xa) 

Index = I 
Do 40 J = I+1,N 
If (X(J) .It. First) then 
First = X(J) 

Index = J 
Endif 
40 continue 

If (Index .ne. I) then 
Tempx = X(I) 

X(I) = X(Index) 

X^ndex) = Tempx 
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Tempy = Y(I) 

Y(I) = Yandex) 
Yandex) = Tempy 
Endif 

30 Continue 
End 


* 


Subroutine Screen 

♦ 

* This subroutine clears the screen, and paints the screen blue. 

* 

Call Init(O) 

Call Clearscreen 
Call Paint(0.0, 0.0,0) 

Call Gscr(l,0,0.0,0.0,0.5) 

Return 

End 

*$include Gkscover.for 
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Appendix D: Examples for the Interpolative Model 


Example 1: Data have been collected for 2090 aluminum in a vacuum environment. The 
data are to be fit to the Forman Equation with closure, fitting only the parameters C and 
n. The constant loading variable is which has a value of 537.6 MPav'mm. The 

frequency of the load is 5 Hz, and the hold time is 0 seconds. The values of Kc and 

AKa, are to be optimized, and the values for the material properties are to be loaded from 
a NASA defined material file. The value of p is to be O.S,and q is to be entered as 1.00 


Note: > indicates keyed input 

Interpolative Model Options 


1) Determine Equation Constants, 
Interpolate Data 

2) Prepare Data File 

3) Quit Program 

Please select an option 


Enter type of units : 

1) U.S. Customary Units (inches, ksi, ksiVin) 

2) SI Units (mm, MPa, MPav^mm) 

>2 


Enter option for Material : 

1) Manual input 

2) Input from NASA/FLAGRO Materials File 

3) Examine data on NASA/FLAGro Materials File 

4) Input from User Defined Materials File (UDMF) 

5) Examine data on User Defined Materials File 
0) Terminate materials entry 

>2 
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Enter one of the codes in brackets [ ] 
to display alloy types or 0 to exit 

NASA/FLAGRO Materials File 


Ferrous Alloys (>=50% Fe) 


Iron (alloy or cast) [A] 
ASTM Spec. Grade Steel [B] 
AISI Steel [C] 
Misc. U.S. Spsc. Grade Steel [D] 
Trade /Common Name Steel [E] 
AISI Type Stainless Steel [F] 
Misc Cres/Heat Resistant Steel [G] 
High Temperature Steel [H] 
Tool Steel [J] 


Non-ferrous Alloys 


1000-9000 Series Aluminum [M] 
Cast and Misc. Aluminum [0] 
Titanium Alloys [P] 
Nickel Alloys/Superalloys [Q] 
Misc. Superalloys [R] 
Copper/Bronze Alloys [S] 
Magnesium Alloys [T] 
Misc. Non-ferrous [U] 


Note: Unless specified properties are for 

1. Welds not heat treated. 

2. 75 deg F Air Environment. 

3. Any Crack Orientation except S-T & S-L. 


>H 


Enter ID Number in brackets to choose a Material Type 
or X to return to main material menu. 


2000 series 


[2] 


5000 series 


[5] 


6000 series 


[6] 


7000 series 


[7] 


>2 




Enter Code in brackets 

to 

choose alloy & HT/Condition . 

2014-T6 A1 

[AD] 

2024-T852 Al 

[EK] 

2014-T651 A1 

[AF] 

2024-T861 Al 

[el] 

2020-T651 A1 

[CB] 

2048-T851 Al 

[FC] 

2024-T3 A1 

[EA] 

2124-T851 Al 

[GC] 

2024-T351 A1 

[EB] 

2219-T62 Al 

[lA] 

2024-T3511 A1 

[EC] 

2219-T851 Al 

[IC] 

2024-T62 A1 

[EG] 

2219-T87 Al 

[IF] 

2024-T81 A1 

[El] 

2324-T39 Al 

[JA] 

2024-T851 A1 

[EJ] 

2090-T8E41 Al-Li 

[PA] 


>PA 


M2PA 2090-T8E41 Al-Li 
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Enter one of the codes in brackets [ ] 

to choose product form, crack ofien. & temperature 

[ 1] M2PA11AB1 Pit & sht; L-T 

>1 

Searching. . . 

Please enter thickness of part (in millimeters) : 

> 10.0 


Material Properties: 


Matl 

No. 

UTS 

YS : Kle 

Klc 

Ak 

Bk 

Thk 

Kc 

1 

586.1 

551.6:1598.4 

1146.7 

1.00 

O 

o 

H 

10.000 

1633.6 


Which equation would you like to use to cuirvefit data : 

1) Forman Equation (with closure) 

2) Forman Equation (without closure) 

3) Paris Equation 

4) Hyperbolic Sine Equation (SINH) 

5) Sigmoidal Equation 

>1 


Would you like to enter the Forman Equation Constants ; 

1) By file 

2) Manually 

>1 


Matl 

No. 


C 


Crack Growth Eqn Constants (closure) 

n : P : <3 : DKo : Rcl :Alpha 



Smax/ 

SIGo 


1 


0.114D-1I: 3.267 : 0 . 50 : 1.00 


4.0 :0.70 :1,50 


0.3 


The following values were determined for Newman's constants : 


AO = 0.3958 

A1 = 0.9255E-01 

A2 * 0.6274 

A3 = -0.1158 


Would you like to use : 

1) Newman Closure (explicit form of Forman Equation) 

2) de Koning/Henriksen equation (Forman variant) 

>2 



Please enter the number of data sets : 

>1 

Please pick one of the following for the data sets: 

1) Constant Stress Ratio 

2) Varying Stress Ratio - Constant Kmax 

3) Varying Stress Ratio - Constant Kmin 

>2 


Load Characteristics for Data Set 1 

Please enter the frequency of the applied load (in Hertz) : 
>5 


Please enter the hold time of the applied load (in seconds) 

>0 

Would you like to : 

1) Optimize the parameters for Data Set 1 

2) Enter the parameters for Data Set 1 

>1 

Determination of Crack Growth Constants 

Press: 1 to enter data from a file, 

2 to enter data by hand. 

>1 


Enter complete file name (maximum of 48 characters) 

>2 09 0TEST.DAT 


Name of data file = 2090TEST.DAT 


Press <return> (or C and <return> to chancre) 

>RTN 

Data exists for Kmax = 537.6 

Press <return> to continue, or C to review and/or edit data 

>RTN 

Save current data on disk? (Y or N) 

Note : New data must be saved now to be used over acrain 
>N ^ 

Curvefit option for Data Set : 1 


n p 

da/dN = C (DKeff) (1 - DKth/DK) 


q 
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[1 - DK/{ (1 - R)Kcr}] 

where DKeff = (1 - So/Smax)/(l - R) 

DKth = DKo (4/PI) arctand - R) 

R = K min / K max 

So & Smax are crack opening & max stresses 
DKo is the fatigue threshold at R=0 

Would you like to curvefit C : 

>Y 

Curvefit option for Data Set ; 1 

n p 

da/dN = C (DKeff) (1 - DKth/DK) 


, q 

[1 - DK/{ (1 - R) Kcr}] 

where DKeff = (1 - So/Smax)/(l - R) 

DKth = DKo (4 /PI) arctand - R) 

R = K min / K max 

So & Smax are crack opening & max stresses 
DKo is the fatigue threshold at R=0 

Would you like to curvefit n : 

>Y 

Curvefit option for Data Set : 1 

n p 

da/dN = C (DKeff) (1 - DKth/DK) 


q 

[1 - DK/{ (1 - R)Kcr}] 

where DKeff = (1 - So/Smax)/(l - R) 

DKth = DKo (4 /PI) arctand - R) 

R = K min / K max 

So & Smax are crack opening & max stresses 
DKo is the fatigue threshold at R=0 

Would you like to curvefit p : 

>N 

Curvefit option for Data Set : 1 

n p 

da/dN = C (DKeff) (1 - DKth/DK) 


q 
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[1 - DK/{ (1 - R)Kcr}] 

where DKeff = (1 - So/Smax)/(l - R) 

DKth = DKo (4/PI) arctand - R) 

R = K min / K max 

So & Smax are crack opening & max stresses 
DKo is t^he fatigue threshold at R=0 

Would you like to curvefit q : 

>N 

Would you like to : 

1) Calculate Kc using NASA Equation 

2) Optimize Kc using Cubic Approximation 

3) Enter Kc Manually 

>2 

Optimum Kc = 0.3297E+09 
p = 0.5000 

Press <return> (or C and <return> to change) 

>RTN 

q = 1.000 

Press <return> (or C and <return> to change) 

>RTN 

Enter 1) Enter AKth manually, 

2) Calculate AKth from DKO 

3) Optimize AKth using Cubic Approximation 

>3 

Optimum AKth = 40.65 

Press <return> to continue 

>RTN 


Alpha = 1.50 

Press <return> (or C and <return> to change) 

>RTN 

Smax/SIGo = 0.30 

Press <return> (or C and <return> to change) 

>RTN 

What confidence interval would you like on the coefficients 
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1) 99.9% 

2) 99.8% 

3) 99.5% 

4) 99.0% 

5) 95.0% 


6) 90.0% 

7) 80.0% 

8) 50.0% 

9) 20.0% 


Please enter selection. 
>6 


Computation of Grack Rate Constants 


S.I. units [mm, MPa, MPaVmm] 

(Crack Closure Considered - Forman Equation) 


Klc « 33.00 Ak = 1.000 BK = 1.000 

YS = 80.00 

Thick = 0.2500 


Kc = 0.3297E+09 


DKo = 56.50 

Sraax/SIGo= 0.3000 Alph^= 1.500 


>RTN 


Press <return> to continue 


Kmax = 
da/dN 

537.6 
Delta K 

0 .2046E-06 

122.8 

0 .4104E-06 

126.1 

0.5594E-06 

129.4 

0.1004E-05 

131.8 

0.8274E-06 

134.4 

0.60221-06 

137.8 

0 .6901E-06 

142 . 0 

0.8955E-06 

146.1 

0.9990E-06 

151.8 

0.2045E-05 

156.5 

0.2255E-05 

161.9 

0.3255E-05 

168.4 

0 .2914E-05 

174.3 

0 .4019E-05 

180.5 


Press <return> to continue 
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>RTN 


da/dN 

Delta K 

0 .4190E-05 

187.1 

0.4424E-05 

195.2 

0.6237E-05 

201.4 

0 . 8037E-05 

207.8 

0 . 1171E-04 

215.2 

0.1286E-04 

222.5 

0.1535E-04 

229.7 

0.1613E-04 

237.6 

0 . 1899E-04 

244 . 8 

0 .2199E-04 

252.8 

0.2506E-04 

260.6 

0 .3016E-04 

267.6 

0 .2826E-04 

275.1 

0.3398E-04 

282.0 

0 .4509E-04 

290.0 

0 . 5364E-04 

303.8 


>RTN 


Press <return> to continue 


Growth Rate Constants: 


C = 

0 .2962E-19 

n = 

6.406 

* p = 

0.5000 

* q = 

1.000 

* indicates a value 

Press <return> to 

continue 

>RTN 


90% Confidence 

Intervals : 

0.6839e-20 s C s 

0.1283e-18 

6 . 120 s n s 

6.692 


input by the user 


Would you like to store the results on a disk ■ 

>N 

The data have successfully been fit to the Forman 
Equation. The coefficient C in the equation has been 
fit to an equation of the form: 

C = A*log ( (1/f) +1) 

where : 
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f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 0.37412839471948E-18 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Equation. The coefficient n in the equation has been 
fit to an equation of the form; 

n = A*log ( (1/f ) + 1) 


where : 

f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 80.899106556360 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Equation. The parameter AKth in the equation has been 
fit to an equation of the form: 

AKth = A*log ((l/f)+l) 

where : 

f = frequency of applied load 

The following values for the. coefficients were 
determined : 

A = 513.38527054904 

Please press <return> to continue 

>RTN 

Comparison of Actual and Predicted Values 


AK 

da/dN (act) 

da/dN (pred) 

% Error 

122.80 

0.2046000E-06 

0.3674280E-06 

79.8 

126.10 

0.4104000E-06 

0.4321204E-06 

5.29 

129.40 

0.5594000E-06 

0.5057272E-06 

9.59 

131.80 

0.1004000E-05 

0.5653662E-06 

43.69 

134.40 

0.8274000E-06 

0.6362456E-06 

23.10 

137.80 

0.6022000E-06 

0.7395572E-06 

22.81 
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142.00 

0.6901000E-06 

0.8853623E-06 

28.29 

146.10 

0 . 8955000E-06 

0.1049106E-05 

17.15 

151.80 

0.9990000E-06 

0.1316127E-05 

31.74 

156.50 

0.2045000E-05 

0.1574851E-05 

22 . 99 

161.90 

0.2255000E-05 

0.1920417E-05 

14 . 84 

168.40 

0.3255000E-05 

0.2413308E-05 

25.86 

174.30 

0.2914000E-05 

0 .2942523E-05 

0 . 98 

180.50 

0.4019000E-05 

0.3592812E-05 

10 . 60 

187.10 

0 .4190000E-05 

0.4403338E-05 

5 . 09 


Please press <return> to continue 


>RTN 





Comparison of Actual 

and Predicted Values 


AK 

da/dN (act) 

da/dN (pred) 

% Error 

195.20 

0.4424000E-05 

0.5584996E-05 

26.24 

201.40 

0.6237000E-05 

0.6644315E-05 

6 . 53 

207.80 

0.8037000E-05 

0.7893201E-05 

1 . 79 

215.20 

0.1171000E-04 

0.9552245E-05 

18.43 

222.50 

0.1286000E-04 

0.1143533E-04 

11.08 

229.70 

0.1535000E-04 

0.1355320E-04 

11.71 

237.60 

0 . 1613000E-04 

0.1619892E-04 

0.43 

244.80 

0.1899000E-04 

0.1892519E-04 

0 . 34 

252.80 

0.2199000E-04 

0.2233087E-04 

1 . 55 

260.60 

0.2506000E-04 

0.2605626E-04 

3 . 98 

267.60 

0 .3016000E-04 

0.2975805E-04 

1.33 

275.10 

0.2826000E-04 

0.3411763E-04 

20 . 73 

282.00 

0.3398000E-04 

0.3850124E-04 

13 .31 

290.00 

0.4509000E-04 

0.4404544E-04 

2.32 

303.80 

0 .5364000E-04 

0.5481080E-04 

2 . 18 


Please press <return> to continue 


>RTN 





Data Set 

Coefficient of Determination 

1 

0.9880 

>RTN 

Would you 

like to interpolate data using the determined 

coefficients 

: 


>N 


Interpolative Model Options 


1) Determine Equation Constants 
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Interpolate Data 


2) Prepare Data File 

3) Quit Program 

Please select an option 


Example 2: Data have been collected for T1-6A1-4V in an air environment. The data are 
to be fit to the Forman Equation without closure, fitting the parameters C, n,p, and q. The 
parameter m is to have a value of 1.0. The constant loading variable is the stress ratio, 
which has a value of 0.1. The frequency of the load is 10 Hz, and the hold time is 0 
seconds. The value of AK* is to be 7.27 ksiv/in, and the values for the material properties 
are to be loaded manually, including the value for Kc of 83.27 ksi\/in. 


Note: > indicates keyed input 


Interpolative Model Options 


1) Determine Equation Constants, 
Interpolate Data 

2) Prepare Data File 

3) Quit Program 

Please select an option 


Enter type of units : 

1) U.S. Customary Units (inches, ksi, ksiVin) 

2) SI Units (mm, MPa, MPaVmm) 

>1 


Enter option for Material : 

1) Manual input 

2) Input from NASA/FLAGro Materials File 

3) Examine data on NASA/FLAGro Materials File 

4) Input from User Defined Materials File (UDMF) 

5) Examine data on User Defined Materials File 
0) Terminate materials entry 

>1 
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Enter Yield Strength: 

> 146.0 

Yield Strength = 146.00 

Enter Plane strain fracture toughness, Klc: 

> 50.0 

Klc = 50.00 

Enter Ak, Bk: 

> 1 . 0 , 0.5 


Ak = 1.00 Bk = 0.50 

Please enter thickness (in inches) : 

> 0.25 

Material Properties: 


Matl 

No. 

YS 

Klc 

Ak 

Bk 

Kc 

1 

146.0 

o 

o 

in 

1.00 

0.50 

62.1 


Which equation would you like to use to curvefit data 

1) Forman Equation (with closure) 

2) Forman Equation (without closure) 

3) Paris Equation 

4) Hyperbolic Sine Equation (SINK) 

5) Sigmoidal Equation 

>2 


m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


(1 


R 


DK/Kcr} 


q 


where : DKth 
R 

DKo 


DKo * (4/PI) * Arctand - R) 

K min / K max 

fatigue threshold @ R=0 


Enter DKo : 
> 3.5 


DKo = 3.50 


Please enter the number of data sets : 

>1 
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Please pick one of the following for the data sets: 

1) Constant Stress Ratio 

2) Varying Stress Ratio - Constant Kmax 

3) Varying Stress Ratio - Constant Kmin 

>1 


Load Characteristics for Data Set 1 

Please enter the frequency of the applied load (in Hertz) : 

>10 


Please enter the hold time of the applied load (in seconds) 
>0 

Would you like to : 

1) Optimize the parameters for Data Set 1 

2) Enter the parameters for Data Set 1 

>1 

Determination of Crack Growth Constants 

Press: 1 to enter data from a file, 

2 to enter data by hand. 

>1 

Enter complete file name (maximum of 48 characters) 

>TO1B.DAT 

Name of data file = TT1B.DAT 

Press <return> (or C and <return> to change) 

>RTN 

Data exists for R = 0.1000 

Press <return> to continue, or C to review and/or edit data. 

>RTN 

Save current data on disk? (Y or N) 

Note : New data must be saved now to be used over again. 

>N 

Curvefit option for Data Set : 1 

m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


,q 

(1 - R - DK/Kcr} 

where: DKth = DKo * (4 /PI) * Arctan(l - R) 
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R 

DKo 


= K min / K max 
: fatigue threshold ® R=0 


Would you like to curvefit C ; 

>Y 

Curvefit option for Data Set : 1 

m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


q 

(1 - R - DK/Kcr} 

where: DKth = DKo * (4 /PI) * Arctan(l - R) 
R = K min / K max 

DKo : fatigue threshold @ R=0 


Would you like to curvefit n : 

>y 

Curvefit option for Data Set : 1 

m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


,q 

(1 - R - DK/Kcr} 

where: DKth = DKo * (4 /PI) * Arctand - R) 
R = K min / K max 

DKo : fatigue threshold @ R=0 


Would you like to curvefit p : 

>Y 

Curvefit option for Data Set : 1 

m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


q 

(1 - R - DK/Kcr) 

where: DKth = DKo * (4 /PI) * Arctand - R) 
R = K min / K max 

DKo : fatigue threshold @ R=0 
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Would you like to curvefit q : 

>Y 

Curvefit option for Data Set ; 1 

m n p 

da/dN = C (1 - R) DK (1 - DKth/DK) 


q 

(1 - R - DK/Kcr} 

where: DKth = DKo * (4 /PI) * Arctand - R) 
R = K min / K max 

DKo : fatigue threshold @ R=0 


Would you like to curvefit m : 

>N 

Would you like to : 

1) Calculate Kc using NASA Equation 

2) Optimize Kc using Cubic Approximation 

3) Enter Kc Manually 

>3 

Enter Kc: 

> 83.27 

Enter 1) Enter Delta Kth manually 

2) Calculate Delta Kth from DKO 

3) Optimize Delta kth Using Cubic Approximation 

>1 

Delta Kth (R = 0.1000 ) = 3.266 

Press <return> (or C and <return> to change) 

>C 

Enter Delta Kth (R = 0.1000 ) : 

> 7.27 

ra = 1 . 000 

Press <return> (or C and <return> to change) 

>RTN 

What confidence interval would you like on the coefficients 


1) 

99.9% 

6) 

90.0% 

2) 

99.8% 

7) 

80.0% 

3) 

99.5% 

8) 

50 . 0% 

4) 

99.0% 

9) 

20.0% 
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5) 95.0% 

Please enter selection. 
>6 


Computation of Grack Rate Constants 
U.S. customary units [inches, ksi, ksiVin] 

(Crack Closure Absent - Forman Equation) 

Klc = 50.00 Ak = 1.000 BK = 0.5000 

YS = 146.0 

Thick - 0.2500 

Kc = 83.27 

DKo = 6.344 


>RTN 


Press <return> to continue 


R = 0.1000 


da/dN 

Delta K 

0.1504E-05 

8.697 

0.1807E-05 

9.049 

0.2390E-05 

9.662 

0.3171E-05 

9.847 

0 .4640E-05 

10.12 

0 .7561E-05 

10.78 

0 .1291E-04 

12.11 

0.1466E-04 

12.41 

0.1372E-04 

13.18 

0.1827E-04 

13.99 

0 .2023E-04 

14.82 

0.2700E-04 

15.63 

0 .3164E-04 

16.39 

0 .3712E-04 

17.39 


>RTN 


Press <return> to continue 


da/dN ; Delta K 
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0.4424E-04 
0 .4232E-04 
0.5009E-04 
0.6002E-04 
0.6819E-04 
0.8099E-04 
0.8054E-04 
0.9050E-04 
0.1006E-03 
0.1130E-03 
0,1785E-03 
0.1681E-03 
0 . 1882E-03 
0 .2449E-03 
0 .2807E-03 
0.2785E-03 
0 .4012E-03 
0.3817E-03 
0.4584E-03 
0,6988E-03 


18.91 

19.36 

21.45 

21.85 

24.46 
25.05 
25.62 
28.07 
28.68 

29.34 
34.31 
34 . 94 
36.11 

38.03 
40,00 

41.85 

44.04 

45.46 

47.34 
51.98 


>RTN 


Press <return> to continue 


Growth Rate Constants ; 


C = 0.2910E-04 

n = 0.2520 

p = 2.257 

q = 2.077 

* m = 1.000 


* indicates a value input by the user 

Press <return> to continue 

>RTN 


90% Confidence Intervals: 

0.7675e-5 s C s O.llOe-3 
-0.2194 s n s 0.7234 
1.961 s p s 2.553 
1.531 s q s 2.623 

Would you like to store the results on a disk ; 

>N 

The data have successfully been fit to the Forman 
Equation. The coefficient C in the equation has been 
fit to an equation of the form: 



C = A*log ( (l/f)+l) 


where : 

f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 0.70299880016317E-03 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Equation. The coefficient n in the equation has been 
fit to an equation of the form: 

n = A*log ( (1/f ) + 1) 

where : 

f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 6.0878585123169 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Ecjuation. The coefficient p in the eqpaation has been 
fit to an equation of the form: 

p = A*log ( (1/f) +1) 

where : 

f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 54.529508388391 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Elation. The coefficient q in the equation has been 
fit to an equation of the form: 

q = A*log ( (1/f) +1) 

where : 
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f = frequency of applied load 


The following values for the coefficients were 
determined ; 

A = 50.178651792115 

Please press <return> to continue 

>RTN 

The data have successfully been fit to the Forman 
Equation. The parameter AKth in the equation has been 
fit to an equation of the form; 

AKth = A*log (d/f)+l) 

where : 

f = frequency of applied load 

The following values for the coefficients were 
determined : 

A = 175.63485719627 

Please press <return> to continue 

>RTN 

Comparison of Actual and Predicted Values 


AK 

da/dN (act) 

da/dN (pred) 

% Error 

8.70 

0.1504000E-05 

0.1228642E-05 

18.31 

9.05 

0.1807000E-05 

0.1887142E-05 

4.44 

9.66 

0.2390000E-05 

0.3291465E-05 

37.72 

9.85 

0.3171000E-05 

0.3771017E-05 

18.92 

10.12 

0.4640000E-05 

0.4520019E-05 

2.59 

10.78 

0.7561000E-05 

0.6509366E-05 

13 . 91 

12.11 

0.1291000E-04 

0.1111932E-04 

13.87 

12.41 

0.1466000E-04 

0.1224719E-04 

16.46 

13.18 

0.1372000E-04 

0.1526259E-04 

11.24 

13.99 

0 . 1827000E-04 

0.1859967E-04 

1.80 

14.82 

0.2023000E-04 

0.2217397E-04 

9.61 

15.63 

0 .2700000E-04 

0.2580105E-04 

4.44 

16.39 

0.3164000E-04 

0.2932364E-04 

7.32 

17.39 

0.3712000E-04 

0.3413289E-04 

8.05 

18.91 

0.4424000E-04 

0.4183280E-04 

5.44 


Please press <return> to continue 

>RTN 


Comparison of Actual and Predicted Values 
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AK da/dN(act) 


19.36 

0 .4232000E-04 

21.45 

0.5009000E-04 

21.85 

0.6002000E-04 

24.46 

0.6819000E-04 

25.05 

0.8099000E-04 

25.62 

0.8054000E-04 

28.07 

0.9050000E-04 

28.68 

0 . 1006000E-03 

29.34 

0.1130000E-03 

34.31 

0.1785000E-03 

34.94 

0.1681000E-03 

36.11 

0.1882000E-03 

38.03 

0 .2449000E-03 

40.00 

0.2807000E-03 

41.85 

0.2785000E-03 


Please press <reti 


>RTN 


da/dN(pred) 

% Error 

0.4420718E-04 

4.46 

0.5585967E-04 

11.52 

0.5821680E-04 

3 . 00 

0.7476238E-04 

9.64 

0.7881556E-04 

2.68 

0.8285450E-04 

2.87 

0.1017637E-03 

12.45 

0.1069109E-03 

6.27 

0.1127023E-03 

0.26 

0.1655161E-03 

7.27 

0.1736471E-03 

3.30 

0.1898270E-03 

0 . 86 

0.2198689E-03 

10.22 

0.2561582E-03 

8.74 

0.2965529E-03 

6.48 


to continue 


Compairison of Actual and Predicted Values 


AK da/dN(act) da/dN(pired) % Error 

44.04 0 .4012000E-03 0 . 3545201E-03 11.64 
45.46 0 . 3817000E-03 0 . 3995588E-03 4.68 
47.34 0 .4584000E-03 0 . 4707966E-03 2.70 
51.98 0 . 6988000E-03 0 . 7325326E-03 4.83 


Please press <return> to continue 

>RTN 


Data Set 
1 

Would you 

>N 


Coefficient of Determination 
0 . 9932 

like to interpolate data using the determined coefficients 


Interpolative Model Options 


1) Determine Equation Constants, 
Interpolate Data 

2) Prepare Data File 

3) Quit Program 

Please select an option 
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Appendix E: Superposition Program Verification - Integration 


To test the accuracy of the linear superposition computer program, it is necessary 
to compare the numerical integration with exact solutions. Exact solutions for each of 
the three major types of loading wave forms (sinusoidal, square, and ramp loading) were 
calculated for a single power law description model of stress corrosion da/dt-K data. 

The integration was calculated over the period of one cycle of loading. As stated earlier: 

=(cK{t)‘'dt (E.l) 

BCC-eX9LCt ^ 

For square wave loading, the program uses the following equations to represent 
the instantaneous stress intensity K(t): 

Kmax (for the percentage of t loaded at maximum, a) 

K(t) = 

K„u„ (for the remaining percentage of r, 1-a) 

Since K in this case is constant for each separate time period, the exact solution to 
Equation E.l is: 

Bcc-exact 


This equation is valid for all C and n. For the tests run, the model was loaded for 50% 
of the cycle, or a = 0.50. 

For ramp loading with a total period, r, the program uses three equations to 
represent K during a load cycle: 
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nr 


(AK/T|o^) * t + K„jn (during the loading period) 

K(t) = (during the hold time) 

- (AK/Tun,o,d) * t (during the unloading) 

Each of these time periods was integrated individually and added to create the final 
solution. Although the solution for the integral during the hold time was the same as for 
the case of the loaded percent of the square wave, the loading and unloading results 
differed for each constant n in the stress corrosion power law equation. For simplicity, 
the exact solutions were calculated and compared with the numerical integration in the 
program for n values of 1,2, 3, and 4. The following are the exact solutions for 
(da/dN),^,^ corrosion Equation E.l for the loading and unloading K(t) functions and the 
given exponent: 


= C' (Kax - 4^) (E.4) 

unload \ ^ I 
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n = 3 
n = 4 

~^jj *” ^ ^unload [‘Lj-2K^=iK*2K^‘&K‘ -K^AK ’ + ] 

^unload V 3 / 

For the test, the model was statically held for 40% of the period, while loading 
and unloading each required 30% of the cycle. 

For sinusoidal loading with period, t, the program uses the following equation to 
represent the loading function: 


Kit) = 


cos + l] + 


Integration of Equation E. 1 for this loading function yielded unique results for each 

exponent calculated: 

n=l 


cfa 

dN 


0 cc-exact 




(E.7) 


(E.8) 


(E.9) 


(E.IO) 


(E.ll) 


(E.12) 
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n=2 


n=3 


da 

scc-exAct 


C X 


Iaki 

I 8 


+ K ^ 

"■avg 


(E.13) 


dN 

gcc-exact 




(E.14) 


n=4 


da 

dN 


= C X 


acc-6xact 


128 


■LK 




zr 2 
^avg 


(E.15) 


For evaluation of the integration in the program, the following variables were set: 

Frequency = 1 Hz 

C constant from Equation III. 5 (arbitrary units) = 1 
Percent reduction of crack growth in unloading =0% 

For all cases examined, precise agreement was observed between the program numerical 

integration and the exact solutions presented inequations E.2 through E.15. Specific 

results are presented in Tables E. 1 through E. 4. Since C was equated to unity, the 

units for and da/dN are meaningless physically. 

Table E.l: Program predictions of (da/dN)^ for n = 1 


V 

1x10-’ 

100,000 

20,000 

100,000 

^min 

0 

0 

10,000 

99,999 

Sinusoidal 

4.9999995x10-® 

49,999.995 

14,999.999 

99,999.50 

Exact Solution 

5x10-® 

50,000 

15000 

99,999.5 

Ramp loading 

7.0x10-® 

70,000.0 

17,000.0 

99,999.70 

Exact Solution 

7 xlO-® 

70,000 

17,000 

99,999.7 

Square loading 

5.0x10^ 

50,000.0 

15,000.0 

99,999.50 

Exact solution 

5x10-® 

50,000 

15,000 

99,999.5 
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Table E.2; Program predictions of (da/dN)^^ conwtai for n = 2 



1x10-* 

100,000 

20,000 

100,000 

^min 

0 

0 

10,000 

99,999 

Sinusoidal 

3.7499996x10" 

3.7499996x10’ 

2.37499986x10* 

9.99990x10’ 

Exact Solution 

3.75x10" 

3.75x10’ 

2.375x10* 

9.9999x10’ 

Ramp loading 

6.0x10" 

6.0 xlO’ 

3.0x10* 

9.999940x10’ 

Exact Solution 

6 xlO " 

6x10’ 

3 xlO* 

9.99994x10’ 

Square loading 

5.0x10" 

5.0 xlO’ 

2.50x10* 

9.99990 xlO’ 

Exact Solution 

5 xlO " 

5 xlO’ 

2.5 xlO* 

9.9999 xlO’ 


Table E.3: Program predictions of (da/dN)^,^, cwrosion for n = 3 


1^ 

1x10^ 

100,000 

20,000 

100,000 

^min 

0 

0 

10,000 

99,999 

Sinusoidal 

3.1249997x10** 

3.1249997x10'^' 

3.9374997x10'" 

9.999850x10'* 

Exact Solution 

3.125x10'* 

3.125x10''' 

3.9375x10'" 

9.99985x10'* 

Ramp loading 

5.50x10'* 

5.50x10'“' 

5.450x10'" 

9.99991x10'* 

Exact Solution 

5.5x10'* 

5.5x10'* 

5.45x10'" 

9.9999099x10'* 

Square loading 

5.0x10'* 

5.0x10'* 

4.50x10*" 

9.99985x10'* 

Exact Solution 

5x10'* 

5x10'* 

4.5x10*" 

9.99985x10'* 


















































































Table E.4: Program predictions of (da/dN)^^ for n = 4 



X 

p 

100,000 

20,000 

100,000 

^min 

0 

0 

10,000 

99,999 

Smusoidal 

2.7343747x10-^* 

2.7343747x10"* 

6.7734369x10*'* 

9.99980x10” 

Exact Solution 

2.734375x10-*** 

2.734375x10*’ 

6.7734375x10*'' 

9.9998x10*’ 

Ramp loading 

5.20x10-*** 

5.20x10*’ 

1.0120x10*’ 

9.999880x10” 

Exact Solution 

5.2x10*** 

5.2x10*’ 

1.012x10*’ 

9.99988x10*’ 

Square loading 

5.0x10*** 

5.0x10*’ 

8.50x10*'* 

9.99980x10” 

Exact Solution 

5x10-^* 

5x10*’ 

8.5x10*'* 

9.9998x10” 
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Appendix F; Examples for the Superposition Model 


This appendix displays the options available to the user in the linear superposition 
program by proceeding through the program as it predicts corrosion fatigue crack growth 
rates. The example is intended to duplicate results produced by Speidel for corrosion 
fatigue in AA 7079-T651 subjected to cyclic loading in a 3.5 % NaCl solution. This 
material-environment combination is well suited for modeling by linear superposition 
because of the great influence of environment on crack growth. Speidel noted this and 
used linear superposition to model the frequency dependence of crack growth^®. 


Crack Growth Modeling 

After an introduction screen is displayed informing the user of the parameters 
which must be input, the next screen allows the user to decide which model the 
co mp uter will use to represent fatigue in an inert environment, (da/dN)fati*ue- Speidel 
used the following equation to model fatigue: 


^ ,6 (M 

cm \ E 


2.1 X 10 "® 


Kc 

Kc-AK 


(F.l) 


Where Kc and AK are in MPav/m. For his calculations, however, he neglected the 
influence of fatigue on crack growth^. Although the current program is not designed to 
use this equation to model fatigue, data points from either a trend line or specific 
da/dN-AK results can be digitized and input to the program. From these data, constants 
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can be determined for any of the modeling equations used by the program (Paris, 
Forman, or Sigmoidal). The Forman Equation with closure was selected for this 
example and proved to be very accurate. The Forman Equation differs from the other 
two models in that the equation constants may alternatively be read from a library of 
material properties on file. It is assumed the user would prefer to read equation 
constants from the library. Once the Forman Equation is selected as the fatigue model, 
the program returns the equation constants from a file. After the constants are specified, 

the program asks if the user would like to change any or all of the constants to values fit 
from input data. 

To begin the search for equation constants, the program asks which units the user 
would like to use. U.S. customery units are AK in ksi/in and da/dN in in/cycle. For S.I. 
units, AK values are in MPa/mm, while da/dN is given as mm/cycle. If all of the equation 
constants are to be fit to data, either selection will have no influence on the output since the 
units of the equation will be determined by the units of the input data. 

On the next screen, shown in Figure F.l, the program displays those materials on file 
for which Forman Equation constants exist. For each class of material listed, further 
specifications of the material will continue until the specific alloy, heat treatment and product 
form are selected. The program will then ask which Forman Equation (closure or non- 
closure) should be used to model data. For this example, the Forman Equation with closure 
was selected. The program then asks for the thickness of the test specimen, for which 10 
mm was input. After displaying material properties, the program will ask if the user would 
like to input the equation constants manually or read them from the library. In this example. 
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the constants were read from the 


library even though they will 
eventually be determined by data points 
from the original fatigue plot. By 
selecting the Forman Equation with 
closure, the user is required to specify 
which type of closure to use. The 
Forman variant was selected. 

The program then asks the user 
whether or not to fit one or more of 
the Forman Equation constants to a set 
of data. After making the appropriate program inputs the user must input the da/dN vs AK 
data digitized from the original graph. After inputing the data and specifying which value 
(K„^, or R) was constant in the data set, the program fits the data to a cubic equation 
of the form: 

— = + AK + C (F.2) 

dN 

where C is a constant. This is shown in Figure F.2. This is done to approximate AK* and 
Kc values for the data set. When the program asks for AK„ and Kc values which represent 
AK* and AKc, respectively, at R = 0, the program displays values from the library file and 
recommended values calculated from the data. Equation F.2 is used to make 
recommendations for the threshold and critical stress intensities from the iiq>ut data. It is 
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suggested that the values recommended 


by the program be used unless the 
operator knows the exact critical 
and/or threshold values to prevent the 
program from having a fatal error. 

After inputing the data, the 
program asks if the data should be 
"flipped". The flipping option allows 
the program to remove some of the 
randonmess of experimental data. The 
random nature of the experimental 
data, when input into the program, can 
result in equations which do not follow 
the general trend of crack growth 
models. Flipping helps the curve 
fitting subroutine to avoid creating 
equations which fit the data but do not 
accurately model crack growth. An 


, . Yod Oivtloii of lh« h th* loKowlfi$ ^u«tk>n: 

<h! . . ■ n n 

; ^ '"r-. 

cTJikck opening 
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2 ) 
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ajKimax 

Pteaee Input your choice^ 

WhaMi your constant vaftM? 

6,0 - ="t 


Figure F.2: Input of data modeled by- 
Forman Equation. 


example of inaccurate modeling is an equation where increasing AK causes decreasing da/dN 
for any AK. Flipping of data should not be done to data sets with few points (< 15) because 
the process deletes points outside of die general trend of the data; the resulting number of 
points after this process may not be sufficient to allow the progr am to fit the Forman 
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Equation to the data. In this situation, 
the data were not flipped. 

Next, the program asks the user 
if each equation constant should be Ht 
to the data or kept at its current value. 
For the example problem, all constants 
were fit to the data, resulting in a 
maximum error of 3.6% between the 
Forman Equation and the input data. 
The program then allows the user to 
graph the input data aixl the fitted 
curve. 

At this point, the program 
begins modeling stress corrosion. For 
simplicity, the power law model was 
chosen. Next, the choice between a 
single and double power law is made, 
followed by the decision to allow the 
program to calculate the constants. 
Stress corrosion data for 7079-T651 in 
3.5% NaCl solution are available from 


the original paper presented by 


ErW t 

ii-iiiliiiiiiii 

Kef- 2t,00 

p^f3CX)00i-<86" 

Eniltr Alphi* 

Alphll ^ T JO 

Smax/^Qo, t«de ef lo Oo^r«*s ^rren% - "0.d4^ 

- OJO 

WouNf you flke ^ fH tor {cufwnHy g.f01000> 

OJOOOOO^ 

Wouidyou^tolKtorf^fcyn«n0y l*00000(^ 

fWicfton Oo^ coofffclwif 
ij 

ib i:oo4$7$da 

MiiiKNurw lrmr» 0.34074 Wiwlimum Porponf Error- 
&lu\dtnS Etovtet)oH’* 0.1742470 

OpfTOleOprv OJWOOO^KWOeO 

Erofi ^'to<uro> to coning* 

'EoftoOfl cwiftofiif 

* p * O.7S08 

*»lndfeeto»pv»lug4Uiv»fitlromrto 



156 


Speidel“. The program was employed to calculate the constants from the input data. The 
program also allows the user to alter the data before importation into the program. One 
method of alteration puts an end of file marker in the data set notifying the computer where 
to stop as it reads in data. Ordering the data from lowest to highest AK values is required to 
insure that the program will not shut down during operation. The conversion of units will 
allow the user to multiply the stress intensity (X) and the crack growth rate (Y) columns by 
separate multipliers. The user must calculate the X and Y multipliers to convert the data 
from current units to the units consistent with inert environment fatigue crack growth rate 
data. 

After the computer imports the stress corrosion cracking data, the program performs 
the flipping routine and calculates the best fit for a double power law relating da/dt to K. 

The two sets of constants are output, as well as the stress intensity where the equation 
changes from using the first set of constants to the second. The program again allows the 
user to view the data graphed with the model. When data are graphed, the program will ask 
if the user would like to write the curvefit data to a file. This option allows the user to 
import the model into a graphing package (Quatro Pro, SigmaPlot, etc.) for later use. This 
option is also available for the fatigue data. 

Output Data Parameters / Calculation 

After modeling the two components of crack growth, the program asks a sequence 
of questions to specify the output data variables for the environmental fatigue crack 
growth rate prediction. These include the constant loading parameter, range of AK 
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values, spacing of AK within this range, and the number of data points. The next screen 
gives the user the opportunity to reduce the calculation of crack growth due to 
unloading. Finally, the type of loading is specified. In this example, sinusoidal loading 
was selected to agree with the loading condition used by Speidel. After inputing the 
frequency, the program predicts (da/dN)„,^ versus AK. The following screen allows the 
user to type in the X and Y units for the output data file, then a location where the 
progr am sends the file. The program then displays the results on the screen, including 
the contribution of stress corrosion to the total fatigue crack growth rate. The program 
then asks if the user would like to run the program again. This will only allow the user 
to alter the loading function; the predefined models and the output parameters will 
remain the same. After selecting no, the program ends. 
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Appendix G: Computer Modeb for Fatigue Analj'sb 


The programs developed by Jonathan Pope, Allen Wilson, and Edward Richey are 
available m executable form, along with the source codes. To tun the programs, place 
the program disk in drive A:\ and type MENU. A tide screen will appear; press 
< RETURN > and a main menu will appear. The main menu allows the user to choose 

between digitizing a set of data, or using either the linear superposition or intetpolative 
models. 

If the user wishes to tun the programs from a drive other than A:\, or from a 

hard drive, the batch file SETWGKS.BAT must be modified. If the user wishes to install 

the programs onto a hard drive, the program files must all be copied into the desired 

directoo- on the hard drive. The batch fiie must then be edited. The file should contain: 
©echo off 

set library = drive-, drive\w!Afor,drive\wgks 
wgksl 

where drive specifies the location of the GKS (Gmphical Kernel System) iibmy files. For 
example, if the user wishes to install the program on the hard drive (disk C:\) in a 

directory called NASA, then drive = C:\NASA, and setwgks.bat would contain: 

©echo off 

set library = C:\NASA;C:\NASA\watfor;C:\NASA\wgks 
wgksl 
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If the user wishes to run tiie program from a floppy disk iii a' drive other than A:\, then 
only the file setwgks.hat must be modified. 

The program files found on the program disk include all of the necessary 
executable files, and the materials files. These files are : 

SETWGKS.BAT - defines the location of die WGKS’ library files 
MENU.EXE - executable form of main menu program, written in QBasic 
NASA.EXE - executable form of iiiterpblative model prb^am 
LINEAR.EXE - executable form of linear superposition model 
DIGITIZE.EXE -executable form of graph digitization program 
ROTUNDA.EXE - executable program which draws main title screen 
WGKS1.EXE - loaJs WAtSOM WgKS intb meihofy 
WATFOR.LIB, WGKS.OB - library files for WGKS 
NASAMFM, NASAMFC, USRMFM, USRMFC - rhatefials files 
The source code listings are also included on the pfogim" disk. The following files are 
needed for each of the executable programs: 

NASA.EXE - Interpolative model 
Nasa.for 
Graph, for 
Filesub.for 
Cfit.for 
Forman, for 
Sig.for 
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Sinh.for 


LINEAR. EXE - linear superposition model 
Linear, for 
Linear2.for 
Inmfil.for 
Graph2.for 
DIGITIZE.FOR 
Digitize, for 
ROTUNDA.EXE 
Rotunda, for 

All of the programs also require the files Gksdefn.for and Gkscover.for. 

The source codes for the computer programs are contained in the compressed file 
SOURCE.EXE located on the program disk. To expand the SOURCE.EXE file, copy 
the file to the directory you wish to use and type SOURCE.EXE. The datafiles listed in 
Appendix A are contained in the compressed file DATAFILE.EXE. Both files are self- 
extracting. 

The executable programs were all complied using the WATFOR-77 compiler 
from Watcom. The compiler only uses conventional memory, thus, in order for the 
program to run, there must be at least 570K of free conventional memory before loading 
WGKS into memory. Urmecessary memory resident programs, such as mouse drivers, 
should be removed. 
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